AFWAL-TR-80-4052 


TIME-DEPENDENT  ENVIRONMENTAL  BEHAVIOR 
OF  GRAPHITE/EPOXY  COMPOSITES 


K  G.  KIBLER 

GENERAL  DYNAMICS  CORPORA  TION 

P.  0.  BOX  748 

FORT  WORTH.  TX  76101 


MAY  1980 


TECHNICAL  REPOl.  r  AFWAL-TR-80-4052 

Final  Report  for  period  1  August  1977  —  31  January  1980 


DTIC 

ELECTE] 

OCT  3  0  1980 


B 


Approved  for  public  release , distribution  unlimited. 


MATERIALS  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related  Government  procurement  operation, 
the  United  States  Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever;  and  the  fact  that  the  government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be  re¬ 
garded  by  implication  or  otherwise  as  in  ary  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture, 
use.  or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (ASD/PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS)  .  At  NTIS,  it  will 
be  available  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


ff  M.  WHITNEY,  Project  g^gineer 
Mechanics  &  Surface  Interactions  Br . 
Nonmetallic  Materials  Uivision 


52 

S.  W.  TSAI, 

Mechanics  &  Surface  Interactions  Br. 
Nonmetallic  Materials  Division 


FOR  THE  COMMANDER 


'F .  D.  CHERRY 
Nonmetallic  Materials  Division 


"If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list,  or 
if  the  addressee  is  no  longer  employed  by  your  organization  please  notify  AFWAL/MLBM, 
W-PA1B,  OH  45433  to  help  us  maintain  a  current  mailing  list". 

Copies  of  tills  report  should  not  be  returned  unless  return  is  required  by  security 
considerations,  contractual  obligations,  or  notice  on  a  specifi:  document. 

AIR  FORC^/56  790/26  5*pt»mt«r  1980  —  200 


4 


CLASSIFICATION  OF  THIS  PACE  (Whan  Dml,  Ent*r*d) 


:iaa»i:f»nwit'.W!«™niai>J«l 


|rNUM0ER 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


RECIPIENT  S  CATALOG  NUMBER 


4.  TITI  F  Qfirf  Sub, I,!.  1  - - - - - - 

f  Time-Dependent  Environmental  Behavior  \ 
of  Graphite/Epoxy  Composites  ~~ 


4.  PERFORMING  ORG.  REPORT  NUMBER 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  4  WORK  UNIT  NUMBERS 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

General  Dynamics 

P.  0.  Box  748  ^  /OV* 

Fort  Worth,  TX  76101  _ 


It.  CONTROLLING  OFFICE  NAME  AND  AOORESS  '*•  REP 

Materials  Laboratory  '  f  jl 

Air  Force  Wright  Aeronautical  LaboratoxLew  u.  number  of  pages 
AFSC ,  Wright-Patterson  AFB ,  OH  45433  |  1^0 _ 


.  MONITORING  AGENCY  NAME  A  AOORESSflf  dltfotoni  from  Controlling  Ottico)  15.  SECURITY  CLASS,  (ot  thin  roport ) 


~2419^j  03  04  (Zj/feX-/ 


l5'l- 


Unclassified 


IS«.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (at  thlt  Rtport) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ot  tho  obmtroct  ontorod  In  Block  20,  If  dlf  for  on  t  from  Roport) 


19.  KEY  WORDS  fConrinu#  on  rooorao  oido  It  nocoaaory  and  tdonttty  by  block  nvmbor) 

Composite,  Viscoelas Licitv ,  Graphite-Epoxy,  Tine-Dependent  Behavior, 
Polymeric  Materials,  Time-Temperature-Humiditv  Superposition, 
Environmental  Creep  Testing. 


ABSTRACT  fConf/nuo  on  tooorao  a/do  if  nocoaaory  and  ‘dontify  by  6/ocJr  mjmbor) 

Experimental  and  analytical  investigations  have  been  conducted  to 
determine  the  individual  and  combined  effects  of  temperature  and 
absorbed  moisture  on  the  time-dependent  mechanical  response  of  uni 
directional  and  cross-plied  graphite-epoxv  composites.  Tensile 
creep  response  was  measured  at  room  temperature  and  four  elevated 
temperatures  (up  to  the  glass  transition  region  of  the  materials) 
for  dry  specimens  and  for  specimens  moisture-saturated  (continued) 


FORM 
I  JAN  73 


EDITION  OP  1  NOV  65  IS  OBSOLETE 


SECURITY  classification  OF  This  page 


SCCuHiTy  CL  ASSt  f\ C  A  T!QN  Qf  ThiS  PACEfH>i*n  Ca/«  ErffcJ) 


l\ 

i 

\j 

at;  75%  and  95%  relative  humidity.  Results  were  obtained  for  two 
material  systems:  ilarmco  T300/5208  and  Hercules  AS/  3502  .  The  stress 
levels  used  were  demonstrated  to  be  within  the  range  of  linear  vis¬ 
coelastic  response.  Superposition  techniques  were  used  to  determine 
compliance  master  curves  iron,  the  results  for  individual  temperature 
and  moisture  condition,  and  quasi-elastic  conversion  was  employed  to 
extract  relaxation  modulus  behavior.  Expansional  strains  due  to 
temperature  and  moisture  were  measured.  From  measured  principal 
properties,  laminate  response  was  predicted  using  Lime-stepped  elas¬ 
tic  analyses  for  the  environments  considered.  Analytical  procedures 
for  transient  and  nonuni  form  temperature  and  moisture  conditions 
were  formulated.  An  analytical  method  for  minimizing  residual 
stresses  during  cure  was  developed,  and  sample  results  for  the  pres¬ 
ent  materials  are  presented.  Observations  and  recommendations  re¬ 
sulting  from,  these  investigations  are  included  in  this  final  report. 


SECURITY  CLASSIFICATION  of  This  PAGE,'W7.an  Oaf*  Filtered) 


FOREWORD 


This  final  technical  report  describes  the  work  performed 
under  Contract  No.  F33615-77-C-5109 ,  Project  No.  2419/0304, 
from  1  August  1977  to  31  January  1980.  The  program  was  spon¬ 
sored  by  the  Materials  Laboratory,  Air  Force  Wright  Aeronaut¬ 
ical  Laboratories,  Air  Force  Systems  Command,  Wright-Patterson 
Air  Force  Base,  OH  45433,  with  Dr.  J.  M.  Whitney,  AFWAL/MLBM, 
as  Project  Engineer. 

Dr.  K.  G.  Kibler,  the  General  Dynamics'  Fort  Worth  Division 
Program  Manager,  acknowledges  the  guidance  of  Dr.  W.  G.  Knauss 
(California  Institute  of  Technology),  and  Drs .  R.  A.  Schaoery 
and  Weitsman  (Texas  A  &  M  University)  in  this  program.  Many 
technical  personnel  at  General  Dynamics  contributed  to  this 
program,  including  Messrs.  T.  S.  Creasy,  J.  E.  Halkias ,  J.  D. 
Reynolds,  W.  S.  Robinson,  J.  H.  Fruit,  G.  W.  Matlock,  R.  L. 
Jones,  E.  D.  Dawson,  W.  S.  Margolis,  J.  W.  Head,  R.  L.  Ratzlaff, 
J.  R.  Eisenmann,  F.  C.  Nordquist,  H.  C.  Hoffman,  and  Ms.  J.  R. 
DeZem. 


m 


:* 


i 


:  ft 


i  j 


i  , 
1  * 
i  ' 


SUCTION 

I 

II 


III 


IV 


TABLE  OF  CONTENTS 


INTRODUCTION 
BASELINE  DATA 

2.1  Materials 

2.2  Specimen  Preparation 

2.3  Environmental  Conditioning 

2.4  Baseline  Tensile  Tests 
TIME- DEPENDENT  RESULTS 

3.1  Experimental  Procedures 

3.2  Uni  iirectional  Properties 

3.2.1  Results  for  (0)^  Coupons 

3.2.2  Results  for  (t45)9g  Coupons 

3.2.3  Results  for  (90) ^5  Coupons 

3.2.4  Summary  of  Shift  Factors 

3.3  Laminate  Properties 

3.4  Residual  Strengths 

3.3  Stress- Strain  vs.  Strain  Rate 
EXPANSIONAL  STRAINS 

4.1  Moisture  Absorption 

4.2  Moisture- Induced  Swelling 

4.3  Thermal  Expansion 

^4  .  *4  ^  UJ I  111  let  L  V  U1  cuu 


V 


PAGE 

1 

6 

6 

6 

10 

12 

19 

21 

2fc 

26 

23 

34 

39 

44 

43 

51 
5  9 

52 
59 
5? 
63 


TABLE  OF  CONTENTS  (Continued) 


SECTION 

PAGE 

V 

ANALYTICAL  METHODOLOGY 

69 

VI 

EFFECTS  OF  NON - 

UNIFORM  MOISTURE  DISTRIBUTION 

i  74 

VII 

CURE  CYCLE  OPTIMIZATION 

80 

VIII 

CONCLUSIONS  AND 

RECOMMENDATIONS 

84 

REFERENCES 

88 

ATTEND IX  A  - 

BASELINE  TENSILE  DATA 

91 

ATTENDIX  B  - 

QUALITY  ASSURANCE  PLAN 

9  8 

ATPENDIX  C  - 

NOTES  ON  VISCOELASTIC 
ANALYSIS  OF  LAMINATES 
by  R.  A.  Schanery 

102 

APTENDIX  D  - 

MINIMIZATION  OF  RESIDUAL 
THERMAL  STRESSES  IN 
CROSS-PLY  LAMINATES 
bv  Y.  Weitsroan 

126 

LI  ST  OF  ILLUSTRATIONS 


FIGURE 


'  PACE 


Pro  ,r;v.  i  Plan 


i  e  s  t  Ma  L  r  1  x 


Time  Dependent  'res  ts  / Predict  ions 

Liquid  Chromatograms  of  Hercules  3502  and  Narmco 
5208 


Coupon  Configurations 
Sample  Stress-Strain  Output 

Retension  of  Matrix-Dominated  Tensile  Properties  at 
Llevated  Temperature  for  T300/5208 


I  UiSOOii  ivti 


:io  for  Unidirectional  (0)^  Coupons 


,  ■  ^ 

’  M 


Superposition  Procedure  -0 

E  9  t  ( t )  at  299°K(RT)  and  339°K(150°F)  for  Dry  T300/5208  22 

-  <L 

Test  Configuration  for  Creep  Experiments  23 

Relaxation  Modulus  for  T300/520S  (FT  Dry),  ( x  ^ )  2  S 

Linear  Viscoelas tici tv  of  T300/5208  (90) 15  :  95%  RH  , 

150°F  ’  27 

Typical  Epp  and  v\2  for  RT  Wet  (0)^  Coupons  29 

Relaxation  Modulus  for  (0)5  AS/3502  29 

Major  Poisson  Ratio  for  (0)g,  AS/  3502  20 

Creep  Compliance  Results  for  (-45) 25  AS/ 3502  30 

Creep  Compliance  Master  Curve  for  Dry  (~t45)  32 

Relaxation  Moduli  fur  Dry  (^45)  32 


LIST  OF  ILLUSTRATIONS  (Continued) 


Poisson  Ratio  for  Dry  (t45)c  T300/5208 

Dry  L22  Master  Curves 

Wet  L22  Master  Curves 

Combined  Wet/Dry  E?2  Master  Curves 

K->2  Master  Curves  for  "Unaged"  Materials 

Long  Terra  Test  Results  Shifted  with  a?  from  t45's 

Dry  Lgg  Result  s  Shifted  with  a-p  from  ±45 's 

Time-Temperature-Humidity  Shift  Factors 

Sx(t)  tor  ( 0/90/145) s  T300/5208 

Me  art-red  6c  Predicted  Compliance  for  Diy  (0/-45)g 
and  (0/90/145) s 

Measured  and  Predicted  Laminate  Response 

Moisture  Induced  Weight  Gain  for  (0)5  Coupons 

Moisture  Induced  Weight  Cain  at  297°K  (RT) ,  754  RH 
for  Various  Laminates 

Moisture  Induced  ’Weight  Cain  at  297°K  (RT)  ,  974  RH 
for  Various  Laminates 

Comparison  of  Moisture  Uptake  in  Different  h-Ply 
Laminates 

Variation  of  "Maximum”  Moisture  Level  with  Relative 
Humidity  at  Pudoit.  Temperature 

Moisture- Induced  Swel ling 

TMA  Schematic  for  Loaded- Co lumn  Thermal  Expansion 
Tests 

Longitudinal  Thermal  Expansion  Coefficient  for 

Both  Materials  Conditioned  at  954  RH 


viii 


LIST  OF  ILLUSTRATIONS  (Continued) 


FIGURE 


PAGF. 


39  . 

Longitudinal  Thermal  Expansion  Coefficients  a  ^ 
for  Dry  Coupons 

c:’ 

40  . 

Transverse  Thermal  Expansion  Coefficients  a >9 
Dry  Coupons 

for 

64 

'  41 . 

Transverse  Thermal  Expansion  Coefficient  a  •>  2 
Coupons  Saturated  at  753)  RF, 

for 

65 

42  . 

Transverse  Thermal  Expansion  Coefficients  a0 1 
Coupons  Conditioned  at  95T  RH 

for 

66 

43. 

Glass  Transition  vs.  Moisture  Content 

6  7 

44  . 

Measured  and  Predicted  Lamina  Properties 

70 

;  45 . 

Exponential  Series  Fit  (*  Decade)  -  t45 

-7  \> 

/  •- 

1  <46  . 

Exponential  Series  Fit  (1  Decade)  -  t45 

7  3 

4  /  . 

i 

i 

Con Jit ioning/Tes t  Sequence  lor  Non-Uniform  Moi 
Content 

s  r  tire 

75 

!  48 . 

Schematic  Procedure  for  Predicting  Mechanical 
with  Non- Uni form  Moisture  Content 

Response 

76 

:  if 

Non-Uniform  Moisture  Tests/Predictions 

7  7 

*1 

1  •  50  ■ 

*  1 

Power  Law  Fits  to  90°  Compliance 

81 

;  r  1 

*  JJ.  . 

Optimized  Cooling  Paths  for  AS/ 3502 

8  2 

7 


4 


LIST  OF  TABLES 


I 

( 


! 


\ , 
i 


i  4 


( . 
i 


TABLE 

I  Typical  Thickness  Variation  in  Test  Coupons 

II  Dry  Resin  Content  and  Specific  Gravity  for 

T300/520S  and  AS/3502  Specimens 

III  Tensile  Strength  Retention  (Relative  to  RT  Dry) 

IV  Residual  Strength  Data  -  AS/ 3502 

V  Residual  Strength  Data  -  T300/520S 

VI  Residual  Strength  Trends 

VII  Conditioning  Parameters 

VIII  Load/Rate  l'arame ters  for  aj/  Measurement s 

IX  Uygro thermal  Constants 

X  RT  Strength  Tests 

9 

XI  Residual  Stresses  (MX/nm) 


x 


•  <*  ,  ^ 

■  "VT 


SECTION  I 


1 .  INTRODUCTION 


The  use  of  graphite- epoxy  composites  for  aerospace  struc¬ 
tural  applications  has  evolved  as  a  result  of  their  mechanical 
properties,  such  as  high  strength- to-weight  ratio,  and  their 
tailorability  of  design,  relative  to  conventional  materials. 

The  composite  constituent  materials  are  usually  assumed  to  re¬ 
spond  to  external  inputs  in  a  linearly  elastic  manner,  and  this 
feature  is  reflected  in  the  common  analytical  methodologies 
used.  It  must  be  recognized,  however,  that  composites  do  ex¬ 
hibit  a  significant  amount  of  time,  temperature,  and  moisture 
dependent  behavior,  particularly  in  service-type  environments. 

In  fact,  the  problem  of  time  influences  the  behavior  of  these 
materials  In  several  ways;  e.g.,  the  time  of  exposure  to  envir¬ 
onment,  the  age  of  the  material,  and  the  time  (duration)  of  a 
particular  loading  condition. 

Although  much  research  has  been  conducted  on  the  visco¬ 
elastic  behavior  of  polymers  and  the  particulate-filled  compos¬ 
ites,  only  limited  investigations  on  the  time- dependence  of 
fibrous  composites  had  been  accomplished  through  1975  (Ref.  1-7, 
e.g.)  Recent  years  have  witnessed  rapidly  expanding  activities 
directed  more  specifically  at  time- dependent  behavior  of  ad¬ 
vanced  composite  materials  systems  of  particular  interest  for 
aerospace  applications  (Ref.  8-11) .  These  wcvks  have  established 
that  although  the  behavior  of  fibrous  composites  is  in  general 
not  linearly  viscoelastic,  there  is  often  a  considerable  range 
of  input  (temperature,  humidity,  load)  where  linear  response 
may  at  least  be  a  very  appropriate  approximation. 

The  overall  objective  of  this  program  was  to  determine 
analytically  and  experimentally  the  coupled  and  uncoupled  effects 
of  temperature  and  absorbed  moisture  on  the  time-dependent  mech¬ 
anical  properties  of  unidirectional  and  laminated  graphite- 
epoxy  composites. 

The  scope  of  these  investigations  included  experimental 
creep  and  stress  relaxation  measurements  on  two  material  systems ; 
the  development  of  analytical  methodologies  to  predict  time- 
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dependent  mechanical  response;  an  assessment  of  the  effects  of 
non-uniform  moisture  distribution  on  mechanical  properties;  and 
the  development  of  an  analytical  procedure  for  optimising  cure 
cycles  to  reduce  residual  stress. 

The  technical  investigations  to  fulfill  these  objectives 
were  conducted  in  six  interrelated  experimental  and  analytical 
tasks,  which  may  be  briefly  described  as  follows: 

o  Specimen  Preparation  -  fabrication,  quality  assurance, 
and  Baseline  Data  environmental  conditioning  of 

unidirectional  and  laminated 
tensile  coupons  of  two  material 
systems;  and  measurement  of 
tensile  s tress  -  strain  curves 
in  the  principal  material  direc¬ 
tions  at  room  temperature  and 
elevated  temperature  for  speci¬ 
mens  dry  and  with  two  prescribed 
moisture  contents. 

o  Time-Dependent  Data  -  measurement  of  creep  and  stress 

relaxation  response  in  the  prin¬ 
cipal  material  direction  at  room 
temperature  and  four  elevated 
temperatures  for  dry  specimens 
and  for  those  conditioned  to  two 
prescribed  moisture  contents, 
followed  by  room  temperature  re¬ 
sidual  strength  tests. 

o  Expansionai  Strains  -  measurement  of  longitudinal  and 

transverse  expansionai  strains 
due  to  temperature  and  moisture, 
including  dry  and  wet  glass  tran¬ 
sition  temperature  determinations. 

o  Analytical  Methodology  -  development  of  analytical  proced¬ 
ures  to  predict  the  time-dependent 
mechanical  behavior  of  laminae  and 
laminates,  including  the  effects 
of  residual  stresses. 

o  Effects  of  Non-Uniforrn  -  experimental  and  analytical  assess- 
Moisture  Distribution  ment  ar.d  demonstration  of  the  ef¬ 

fects  of  non  uniform  moisture  dis¬ 
tribution  . 


o  Cure  Cycle  Optimization  -  development  of  an  analytical 

procedure  to  minimize  residual 
stresses  induced  during  the 
cure  cycle  . 

These  investigations,  and  their  inter-relationships,  are 
schematically  outlined  in  Figure  1. 

The  present  work  summarizes  these  investigations  to  obtain 
generic  information  on  the  time-dependent  mechanical  behavior 
of  graphite-epoxy  composites  used  in  aerospace  structures.  Re¬ 
sults  are  presented  for  two  material  systems:  Karmco  T300/5208 
and  Hercules  AS/3502,  both  presently  used  in  F-16  aircraft  pro¬ 
duction.  Experimental  and  analytical  results  are  described  for 
the  individual  and  coupled  effects  of  temperature  and  absorbed 
moisture  on  the  tensile  static  and  creep  properties  of  unidirec¬ 
tional  and  cross-plied  composites,  Mechanical  response  was  mea- 
sureu  at  room  temperature  and  several  elevated  temperatures  (up 
to  the  glass  transition  region  of  the  materials)  for  dry  specimens 
and  for  specimens  saturated  at  75%  and  95%  relative  humidity. 

Superposition  techniques  were  used  to  determine  time-tem¬ 
perature-moisture  master  curves  from  the  results  for  individual 
temperature  and  moisture  conditions.  These  experimental  char¬ 
acterizations  were  then  used  with  conventional  or  slightly  modi¬ 
fied  analytical  techniques  to  predict  time- dependent  laminate 
response,  in  order  to  ascertain  the  magnitude  and  importance  of 
time,  temperature,  and  humidity  effects  on  structural  composites. 

Figure  2  briefly  summarizes  the  test  environments  used  for 
static  and  creep  testing,  while  Figure  3  demonstrates  the  over¬ 
all  features  of  the  test/prediction  scheme  employed. 

In  the  following  sections  are  described  the  procedures  and 
results  for  each  of  the  tasks  indicated  previously  in  Figure  1. 


OBJECTIVE:  Determine  analytically  and  experimentally  the  coupled  and  uncoupled  effects 
of  temperature  and  absorbed  moisture  on  the  time -dependent  properties  of 
g,rat  hite-epoxy  ccnposites. 


.Analytical  Procedure  to 
Minimize  Residual  Stress 


Figure  2  Test  Matrix 


(T30U/5208  &  AS/3502) 


Figure  3  Time-Deperident  Tests/Predictrons 


SECTION 


II 


2 .  BASELINE  DATA 

Summary  details  are  provided  in  this  section  of  the  mater¬ 
ial  characteristics;  specimen  preparation,  inspection,  and  con¬ 
ditioning;  and  the  baseline  tensile  s tress - s train  results.  Com¬ 
plete  baseline  tensile  data  are  rep^  ted  in  Appendix  A  to  this 
report . 


2.1  Materials 

All  experiments  performed  in  this  program  utilized  coupons 
fabricated  from  two  material  systems:  T300/5208  (Narmco)  and 
AS/3502  (Hercules).  Both  materials  are  currently  used  by  General 
Dynamics  in  F-16  production.  As  will  be  seen  from  the  results  to 
be  presented  throughout  this  report,  these  materials  systems  are 
similar  in  their  physical  characteristics  and  mechanical  behavior 

Figure  A  shows  a  typical  liquid  chromatographic  scan  of  the 
two  resin  systems  and  demonstrates  their  similar  chemical  char¬ 
acteristics.  The  as-received  Hercules  material  tends  to  exhibit 
less  of  the  reaction  product  than  Narmco ' s  material.  This  sug¬ 
gests  thac  Hercul.es  may  use  a  somewhat  greater  amount  of  the 
high-viscosity  minor  epoxy  and  thus  does  not  need  to  advance  the 
material  as  far  as  Narmco  does  to  achieve  the  same  viscosity. 


2.2  Specimen  Preparation 

All  specimens  and  coupons  used  in  this  program  were  accep¬ 
ted,  processed,  and  inspected  in  accordance  with  standard  General 
Dynamics'  material  and  processing  specification  as  briefly  out¬ 
lined  in  Appendix  B.  Except  for  thermal-  and  moisture-expansion 
specimens  and  moisture  travelers,  the  coupons  were  of  the  de¬ 
signs  shown  in  Figure  5  All  coupons  mechanically  tested  were 
instrumented  with  2  axial  gages  on  opposite  sides  of  the  speci¬ 
men  as  indicated  in  Figure  5.  In  each  set  of  replicates  for 
(0)6  and  (tA  5)c  tests,  a  transverse  gage  was  also  used  for  ob- 
taininc  Poisson  ratio  data. 
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Figure  4  Liquid  Chromatograms  of  Hercules  3502  and  Narmco  5208 


Figure  5.  Coupon  Configurations 


A  summary  of  the  test  coupons  for  mechanical  response  and 
their  utilization  is  indicated  below: 


o  Unidirectional  Composite  Properties 


-  Baseline 

-  Creep/Relaxation 

-  Residual  Strength 

o  Laminated  Composite 

-  Baseline 

-  Creep/ Relaxation 

-  Residual  Strength 


.  (0)6 

(90)15 
(±45) 2s 

Properties 

^  (0/±45)s 

(90/±45)s 

J  (0/90/±45)s 


o  "Thick"  Specimen  Properties 


-  Effects  of  Lon-Uniform 
Moisture  Eistrition 


> 


(90/+45), ~ 
<±«>6S 


All  materials  were  cured  per  General  Dynamics'  specifica¬ 
tion  which  culminates  with  a  2-hour  hold  at  453°K  (355°F)  with 
5.9  x  105  N/m^  (85  psi)  pressure,  followed  by  a  subsequent  post¬ 
cure  for  6  hours  at  453°K. 


As  previously  noted,  the  quality  assurance  provisions  for 
materials  and  processing  were  conducted  as  outlined  in  Appendix 
B.  Inspection  of  specimens  per  that  plan  did  reveal  irregular¬ 
ities  in  certain  sets  of  (9C/^45)c  and  0°  coupons .  New  coupons 
were  subsequently  fabricated.  The  tensile  data  presented  in 
Appendix  A  indicate  that  problems  still  existed  with  some  of 
the  unidirectional  coupons .  Since  the  0°  coupons  are  highly 
sensitive  to  end  tab  misalignments,  we  believe  that  end  grip 
features  account  for  most  of  the  test  variabilities.  Since 
the  primary  purpose  of  the  baseline  tensile  data  was  to  aid  in 
chocsing  s tress / s train  levels  for  the  time-dependent  tests,  the 
end  grip  problems  were  likely  not  a  factor  for  the  remainder  of 
the  program  since  subsequent  tests  we/e  conducted  well  below 
the  region  of  tensile  ultimate  properties. 

Tabic  I  summarizes  the  typical  thickness  variations  noted 
for  four  differnt  layups  of  each  material.  Our  past  experience 
had  indicated  that  more  thickness  variation  occurs  in  systems 
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using  tihc  AS  fiber  compared  to  T300,  and  this  observation  was 
borne  out  by  the  measurements  given  in  the  table.  Both  mater¬ 
ials  were,  however  ,  withir.  our  nominal  ply  thickness  requirement, 
0.14  +  .001  cm  (.0055  t  .0004  in.). 

TABLL  1.  Typical  Thickness  Variation  In  Test  Coupons 


Table  II  summar ices  resin  content,  specific  gravity  and 
fiber  volume,  Vf,  for  the  cured  coupons.  In  all  cases  except 
those  labeled  "(O)^  large",  the  measurements  were  conducted  on 
test  tabs  cut  from  the  panels  used  for  tensile  coupon  fabrica¬ 
tion.  The  large  (0)£  panels  were  especially  fabricated  to  pro¬ 
vide  the  sizable  number  of  unidirectional  specimens  (normally 
3  cm.  square)  used  for  moisture-induced  swelling  measurements. 
Two  separate  determinations  of  resin  content  were  made  on  24- 
ply  tabs  to  obtain  an  estimate  of  the  through- the- thickness 
resin  content  gradient. 


2.3  Environmental  Conditioning 

For  each  coupon  type  tested  in  this  program,  results  were 
obtained  under  dry  and  wet  (saturated  to  uniform  moisture  con¬ 
tent  at  two  different  humidities)  conditions.  Coupons  for  dry 
tests  were  initally  desiccated  for  several  days  at  moderately 
elevated  temperatures  (322-339°K),  tr.en  stored  in  desiccant  at 
RT  un  til  c  e  s  t . 
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TABLE  II.  DRY  RESIN  CONTENT  AND  SPECIFIC  GRAVITY  FOR  T30C/5208  AND  AS/  3502  SPECIMENS 
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Wet  coupons  were  moisture-conditioned  to  saturation  at 
339°K  (150°F)  in  75?0  and  95  RH .  The  latter  humidity  was 
selected  to  provide  the  greatest  amplification  of  moisture 
effects  without  immersion  conditioning.  The  75%  RH  represents 
a  worst-case  service  condition.  These  specimens  were  supported 
on  racks  in  the  humidified  stagnant  air  space  above  saturated 
salt  solutions  of  lead  nitrate  (  95%  RH)  and  sodium  chloride 
(  75%  RH) .  Coupon  end  tabs  were  wrapped  with  aluminized  tape 
to  retard  moisture  ingress  during  conditioning.  As  previously 
indicated  in  Figure  2,  typical  saturation  moisture  contents 
for  these  two  conditions  were  1 . 0 - 1 . 2%  and  1.4-1. 5%  by  weight. 
Section  IV  details  the  moisture  absorption  characteris tics  of 
these  materials. 


2.4  Baseline  Tensile  Tests 


After  appropriate  environmental  conditioning  was  completed, 
coupons  to  be  tested  were  strain-gaged,  then  reconditioned  (if 
wet)  as  necessary.  Hicro-measurettients  CEA-06- 250UW- 350  foil 
gages  (350  ohms),  .64  cm  (.25  in)  long,  were  bonded  to  the  cou¬ 
pons  with  Eastman  910  adhesive  for  RT  tests,  and  wich  M-Bond  600 
fo1'  all  elevated  temperature  and  wet  tests.  This  latter  adhe¬ 
sive  requires  buying  at  450  i\  (350  F)  for  %  hour  f 


cure  . 


ti¬ 


er 


this  bake,  the  wet  coupons  were  then  reconditioned  for  three  days 
in  the  appropriate  humidity  chamber  to  re-establish  the  outer- 
plys’  moisture  level  lost  in  the  bake. 


All  baseline  tensile  tests  were  performed  with  an  Instron 
tester  controlled  by  a  Hewlett-Packard  Model  B050B  computer. 

The  system  includes  a  10-channel  scanner  for  strain  gage  reading, 
a  10-channel  actuator  scanner  of  control  of  the  load  machine, 
and  a  B  6c  F  signal  conditioner  for  the  strain  gage  inputs.  The 
load  and  strain  inputs  are  processed  by  the  computer  and  output 
via  thermal  printer  and  plotter.  A  typical  sample  output  for 
this  program  is  indicated  in  Figure  6. 


Due  to  the  relatively  short  time  required  for  typical  ten¬ 
sile  tests  (several  minutes),  provision  was  not  made  to  retard 
moisture  egress  from  the  wet  specimens  during  these  tests.  In 
examining  the  wet  baseline  tensile  data,  then,  one  must  note 
that  some  coupon  dryout  does  occur  during  the  test. 


For  each  moisture  condition  (dry  or  wet),  baseline  tensile 
tests  were  conducted  at  room  temperature  (RT)  and  one  elevated 
temperature  (ET) .  The  elevated  temperature  was  chosen  to  be 
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near  the  (dry  or  wet)  glass  transition  temperature  (Tg)  of  the 
materials,  and  to  cause  a  decrease  of  matrix-dominated  tensile 
properties  to  about  70%  of  the  RT  values .  Figure  7  demonstrates 
that  this  goal  was  achieved  at  our  selected  temperatures  of 
A36°K  (325bF)  dry  and  372°K  (210°F)  wet. 

A  tabular  summary  of  tensile  strength  retention  relative 
to  RT  dry  conditions  is  given  in  Table  Ill.  One  notes  that 
slight  problems  still  persisted  with  even  the  refabricated  sets 
of  (0)  and  (90/ 1  <4 5 )  coupons,  since  the  RT  dry  values  'were  low 
compared  to  other  test  conditions.  In  general,  however,  the 
reduction  in  ultimates  due  to  temperature  and  moisture  content 
are  fairly  typical.  Subsequent  Lime-dependent  tests  (Section 
III)  used  stress/strain  levels  which  were  typically  20-30%  of 
ul timatcs . 


A  complete  listing  of  the  baseline  tensile  properties  for 
all  coupons  for  this  task  and  the  following  one  are  presented 
in  Appendix  A. 


Two  results  related  to  time  dependence  were  obtained  during 
baseline  testing,  and  we  will  comment  briefly  on  those  findings. 
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all  the  baseline  tensile  tests.  These  plots  of  transverse  strain 
vs.  longitudinal  strain  demonstrate  that  v -p  is  essentially  con¬ 
stant  over  the  entire  range  of  temperature  and  humidity  (and 
strain)  used  here.  This  provides  strong  evidence  early  in  the 
program,  for  the  time- independence  of 


During  the  later  stages  of  baseline  tensile  testing  of  wet 
coupons,  we  also  initiated  a  modification  of  the  tests  to  obtain 
stress -relaxation  data.  In  these  cases,  we  programmed  the  In- 
stron  tester  to  ramp  at  constant  crosshead  speed  to  a  preselected 
strain  level  (typically  20%  of  ultimate)  and  hold  there  for  30-60 
minutes.  During  the  hold,  all  stress/strain  values  were  printed 
at  regular  intervals  to  observe  the  relaxation  of  the  coupons. 
After  the  selected  hold  time  was  complete,  the  ramp  to  failure 
was  completed. 


Two  interesting  results  were  obtained  from  these  tests: 

(1)  The  hold  period  appeared  to  have  no  influence  on  the  even¬ 
tual  ultimate  properties  measured,  so  no  baseline  data  was  lost 
in  this  process; and  (2)  we  were  able  to  obtain  stress-relaxation 
results  for  the  unidirectional  coupons  which  were  adequate  to 
verify  the  t  ime-  independence  (as  expected)  of  and  V-^n .  On 

Liit  bus  is  of  these  xesui  ls  ,  then ,  wu  uiu  no  l  perform  anv  creep 
compliance  measurements  on  the  (0) ^  coupons.  The  Ep^(t)~E^^ 
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Figure  8.  Poisson  Ratio  for  Unidirectional  (0)g Coupons 


results  from  these  hold  experiments  are  presented  in  Section  III. 
One  should  note  however,  that  those  tests  on  wet  coupons  were 
conducted  at  RT  with  no  provision  for  moisture  maintenance  in 
the  coupons.  Even  though  the  measured  moisture  loss  in  these 
thin  specimens  was  about  20%,  the  E-q  and  principal  moduli 

appear  constant. 

One  word  of  caution  should  be  interjected  regarding  the  use 
of  such  stress-relaxation  data.  Proper  consideration  must  be 
given  to  the  relaxation  effects  which  occur  during  the  irrtial 
ramp  to  the  constant  (hold)  strain  value.  This  is  readily  accom¬ 
plished  by  not  including  the  data  obtained  during  the  early  stages 
cf  the  hold  interval.  That  is,  the  relaxation  which  occurs  during 
the  first  time  interval,  approximately  three  times  the  duration  of 
the  ramp,  is  not  representative  of  the  true  stress  relaxation. 


SECTION  111 


3.  TIME-DEPENDENT  RESULTS 

The  objective  of  this  task  was  to  measure  the  time-depend¬ 
ent  mechanical  response  for  both  material  systems  under  a  var¬ 
iety  of  temperature/humidi ty  conditions.  As  previously  indica¬ 
ted  in  Section  I,  these  tests  were  conducted  on  dry  coupons  and 
on  coupons  conditioned  to  a  uniform  moisture  content  at  73%  and 
95%  RH .  At  each  of  these  humidity  conditions,  five  test  tem¬ 
peratures  were  used,  from  room  temperature  up  to  the  glass  tran¬ 
sition  region  for  the  materials. 

Although  most  of  the  experimental  results  in  this  section 
are  presented  as  time- dependent  moduli,  the  majority  of  our  ex¬ 
periments  were  constant  load  creep  compliance  measurements. 
Since  we  set  out  to  obtain  a  linearly  viscoelastic  charac teriza- 
tion  of  the  materials,  we  indeed  found  that  c uas i - e 1  as t ic  conver¬ 
sion  (i.e.,  E(t)  ^  S(t)“l)  gave  very  good  results  as  expected, 
so  the  moduli  presented  here  are  reciprocal  compliances. 

From  these  creep  compliance  results  for  the  various  layups 
tested,  we  thus  obtained  both  unidirectional  (principal)  proper¬ 
ties  and  laminate  properties.  En  ^  and  1^2  were  obtained  from 
tests  on  (0)^  coupons;  E22(t)  from  (90)^5  coupons;  and  G-^Ct) 
from  Eft(t)  and  Vy.,,  ,  the  longitudinal  modulus  and  Poisson's 

ratio,  respectively,  of  a  (^45)oc;  coupon. 

Front  the  individual  creep  measurements  on  each  layup  at 
various  temperature/hunh dity  conditions,  master  curves  were 
constructed  using  the  procedure  diagranuned  in  Figure  9  .  The 
individual  temperature/humidi ty  results  were  horizontally  shif¬ 
ted  on  a  log  compliance-log  time  plot  until  a  smooth  "master 
curve  was  obtained.  Since  our  tests  were  all  conducted  in  lin¬ 
ear  regions  of  material  response,  and  since  this  procedure  gen¬ 
erally  yielded  reasonable  and  repeatable  results,  we  made  no 
attempt  at  vertical  shifting  of  the  data.  Although  others  (Ref. 

8  )  have  performed  vertical  shifts  based  on  arguments,  e.g., 
from  rubber  elasticity,  we  saw  no  compelling  reason  to  do  so  in 
this  case,  since  our  data  were  quite  tractable  without  such  shift 
ing. 
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In  the  following  subsections,  we  provide  details  about  our 
experimental  procedures  and  present  specific  experimental  results 
for  the  coupons  tested.  Some  of  these  results  do  relate  closely 
to  our  Analytical  Methodology  task.  We  reserve  for  that  section 
(V)  tiiose  pecifics  which  treat  the  prediction  of  time-dependent 
response: . 

3.1  Experimental  Procedures 

Tests,  priniariLy  tensile  creep,  were  conducted  on  all  coupons 
under  the  temperature  and  moisture  conditions  previo1  sly  dia¬ 
grammed  in  Figure  2-  As  indicated  by  the  dotted  square  in  that 
figure,  we  initially  planned  to  use  4G4°K  (375°F)  as  our  highest 
dry  test  temperature.  After  conducting  some  early  testing,  at 
that  temperature  and  experiencing  various  premature  coupon  fail¬ 
ures,  we  lowered  that  selected  test  temperature  to  436°K  (325°F), 
somewhat  below  the  dry  Tg  of  the  material. 

Early  in  the  program  we  experimented  with  various  types  of 
time-dependent  testing  in  order  to  compare  the  respective  results. 
Figure  10  shows  a  typical  example  of  such  results.  Using  three 
different  90°  coupons  of  T300/5203,  we  measured  both  stress  re¬ 
laxation  (by  two  different  methods)  anc  creep  compliance.  These 
data,  all  plotted  as  moduli,  are  compared  in  Figure  10.  Tor  the 
tests  shown,  the  largest  disagreement  in  results  was  about  4e; . 

The  choice  of  experimental  method  and  apparatus  for  the  bulk 
of  the  time-aependeni  tests  was  driven  by  several  requirements: 
total  number  of  specimens  to  be  tested,  number  of  replicates  per 
test,  and  maintenance  of  moisture  in  the  wet  specimens  during 
test.  The  experimental  creep  apparatus  i'  shown  schematically 
in  Figure  11. 

This  test  set  up  utilized  a  whif fie- tree  fixture  to  accom¬ 
plish  simultaneous  loading  of  four  specimens  (3  replicates  plus 
one  dummy) .  The  strain  in  each  of  the  replicate  specimens  was 
monitored  with  two  strain  gages  on  opposite  sides  of  the  speci¬ 
men.  The  average  of  the  two  gage  outputs  for  each  specimen  was 
plotted  on  a  strip  chart  recorder.  For  analysis,  the  plotted 
strains  from  the  three  specimens  were  averaged,  and  the  average 
tabulated  vs.  time.  In  the  case  of  t45  cou,.j;\= ,  one  of  the  rep¬ 
licate  specimens  was  some Lines  monitored  with  an  additional  gage 
(transverse)  in  order  to  extract  the  Poisson  ratio  Vxv- 


><t)  at  299°K  <RT)  and  339°K  (150°F)  tor  Dry  T300/5208 


For  dry  tests,  an  oven  was  placed  about  the  fixture,  with 
a  tray  of  desiccant,  and  the  temperature  was  continuously  moni¬ 
tored  by  a  thermocouple  placed  on  one  of  the  specimens . 

For  wet  tests,  we  accomplished  a  somewhat  innovative  pro¬ 
cedure  whereby  the  specimens  were  continuously  immersed  in  liquid. 
Immersion  heaters  and  impellers  were  used  for  temperature  uniform¬ 
ity.  One  of  the  paramount  problems  existing  in  usual  wet  test 
techniques  is  the  maintenance  of  moisture  in  previously  condi¬ 
tioned  wet  specimens  and/or  accounting  for  dryout  during  the  test. 
This  is  not  an  insurmountable  problem  in  the  case  of  tensile  tests, 
where  the  test  times  are  typically  of  a  few  minutes  duration.  In 
the  present  case  however,  our  planned  tests  were  to  be  100  minutes, 
so  moisture  egress  was  a  major  consideration. 

As  discussed  in  Section  IV,  the  wet  glass  transition  temper¬ 
atures  for  these  coupons  were  all  measured  to  be  in  the  vicinity 
of  the  boiling  point  of  water.  This  fact  allowed  us  to  immerse 
the  coupons  in  water  mixtures  and  still  test  at  temperatures  very 
near  the  glass  transition  range.  The  immersion  fluids  selected 
for  these  tests  were  mixtures  of  ethylene  glycol  and  water.  The 
mixtures  were  selected  on  the  basis  of  Raoult's  law  (Ref.  12,  e.g.) 
tc  provide  effective  vapor  pressures  of  water  equivalent  to  that 
at  which  the  coupons  were  initially  conditioned,  i.e.,  75%  or  95% 

RK .  Such  a  procedure  may  be  expected  to  work  reasonably  well  when 
the  liquid  added  to  the  water  has  a  high  boiling  point. 

During  all  the  creep  testing,  the  moisture  content  of  the 
specimens  was  monitored  by  periodic  weight  measurement  of  a  nois- 
Luie  traveler  specimen,  of  the  same  type  as,  and  submerged  with 
the  creep  coupons.  During  a  typical  test  duration  on  one  set  of 
replicates,  three  days,  the  traveler  indicated  only  slight  moisture 
pickup  (uMMlly  less  than  10%  increase). 

Strain  gages  used  during  these  immersion  tests  were  the  same 
type  as  described  previously  in  Section  11.  For  these  immersion 
tests,  however,  we  coated  the  gages  with  Aqua  Seal,  an  elastomeric 
product  sold  by  Kearney  National  ,  Inc.  Although!  this  does  not 
inhibit  moisture  uptake,  it  diu  retard  moisture  effects  on  the 
gages  for  a  sufficient  time  to  cor.plere  a  given  veries  of  tests. 

In  separate  tests  conducted  on  aluminum  coupons  gaged,  coated,  and 
immersed  in  the  same  manner  as  our  graphite-epoxv  specimens ,  we 
verified  that  reliable  strain  readings  were  being  obtained.  Al¬ 
though  gage  problems  would  typically  initiate  after  5-10  days  of 
testing,  this  was  well  outside  our  usual  3-day  test  sequence  for 
five  temperatures  on  one  bee  of  bpeeii.ieus . 


The  typical  creep  testing  sequence  consisted  of  100  minutes 
of  creep  at.  a  given  temperature,  followed  by  a  recovery  period 
with  no-load  while  the  temperature  was  raised  and  stabilized  to 
the  next  value.  There  are  two  time  intervals  during  such  ex¬ 
periments  which  should  be  considered:  the  time  immediately  fol¬ 
lowing  load  application,  and  the  recovery  time  between  success!  e 
loadings . 

In  the  first  case,  one  should  note  that  true  compliance  mea¬ 
surements  are  obtained  only  after  a  time  interval , typically  3-5 
times  the  loading  interval,  due  to  relaxation  during  loading. 

Our  initial  creep  strain  readings  were  recorded  at  0.5  minutes 
after  loading.  The  100  minute  creep  duration  thus  gives  us 
slightly  over  two  decades  of  response,  which  was  acceptable  for 
later  forming  master  curves  since  we  obtained  data  at  five  separ¬ 
ate  test  temperatures. 

In  the  second  case,  the  recovery  time  should  normally  be 
about  three  times  the  creep  interval  for  the  specimen  to  approach 
its  initial  condition.  This  feature  is  readily  ascertained,  how¬ 
ever,  from  successive  identical  creep  experiments.  For  our  pur¬ 
poses,  we  allowed  the  coupons  to  recover  for  a  time  equal  to  the 
creep  interval  (at  the  creep  temperature')  then  raised  the  tempera¬ 
ture  to  tne  next  higher  test  value.  by  the  f'me  the  chamber  had 
equilibrated  at  the  higher  temperature ,  (e.g.,  60-100  minutes), 
the  coupons  had  sufficiently  recovered  for  the  next  test. 

Since  our  intent  for  this  program  was  to  obtain  a  complete 
linearly  viscoelastic  characterization  for  these  materials,  we 
devoted  considerable  attention  to  determining  the  cycle- to- cycle 
repeatability  of  the  experiments  and  the  range  of  linear  behavior. 
During  this  process  we  assessed  the  possible  requirements  for 
mechanical  "break-in"  of  the  coupons  prior  to  actual  creep  exper¬ 
iments  . 

Various  workers  (e.g..  Ref.  3  and  9  )  have  noted  the  need 
to  condition  these  fibrous  composites  mechanically  before  attempt¬ 
ing  to  obtain  reproducible  time-  dependent  data.  V.'e  examined  such 
requirements  by  low  rate  strain  cycling  (  0.057,  per  second)  hys¬ 
teresis ,  and  by  cycle- to-cycle  (constant  load;  creep  recovery  re¬ 
sults.  Although  in  some  cases  a  different  response  was  noted  for 
the  first,  cycle  compared  to  subsequent  cycles,  this  was  not  uni¬ 
formly  the  case.  We  adopted  a  precautionary  procedure,  therefore 
in  which,  prior  lc  the  creep  experiments,  the  coupons  were  sub¬ 
jected  to  3-5  load/unload  cycles  H  min  each)  using  a  107c  higher 
load  than  planned  for  the  actual  creep  experiment.  Following 
this  cycling,  a  30-minute  recovery  period  was  allowed  piiui  to 
initiation  of  the  creep  experiments. 
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A  typical  cycle- to- cycle  repeatability  test  is  indicated  in 
Figure  12,  where  we  also  assessed  the  range  over  which  the  tensile 
relaxation  modulus  of  a  +45  coupon  remains  constant.  For  a  strain 
range  of  0. 1-0.2%,  the  variation  between  maximum  and  minimum  Ep  at 
a  given  time  is  about  5%.  Such  a  range  is  probably  within  the 
typical  specimen  to  specimen  variability.  If  we  include  the  Ep 
data  at  0.3%  strain,  the  variation  between  minimum  and  maximum 
results  increases  only  slightly  to  about  7%.  If,  however,  we 
examine  the  ER(t)  results  over  a  strain  range  of  0.1-0. 4%,  the 
variation  at  a  given  time  increases  to  around  15%,.  These  results 
would  tend  to  indicate  that  for  RT  dry  specimens,  the  material  is 
linearly  viscoelastic  up  to  about  0.3%  strain,  which  is  roughly 
30, s  of  the  strain  ultimate  under  these  conditions.  The  stress 
corresponding  to  this  strain  level  is  about  60  MN/m^  (8.75  ksi)  , 
which  corresponds  to  about  35-40%  of  ultimate  stress. 

Stress  levels  for  all  creep  tests  were  set  at  typically  25-  35!, 
of  ultimate  values  for  each  temperature/moisture  condition.  For 
a  given  sequence  of  test  temperatures  at  a  particular  moisture  con¬ 
dition,  the  creep  load  was  systematically  lowered  at  each  higher 
test  temperature  in  order  that  the  creep  strains  remained  approx¬ 
imately  constant. 

Verification  of  linear  response  was  obtained  by  superposition 
of  creep-recovery  results  (Ref.  7  )  as  shown  for  a  typical  case  in 
Figure  13  where  a  comparison  is  made  between  measured  and  extrapol¬ 
ated  compliance  at  equivalent  times  in  the  respective  creep  and 
recovery  periods . 


3.2  Unidirectional  Properties 

In  this  subsection,  we  describe  the  experimental  results  from 
which  principal  properties  for  these  two  material  systems  are  de¬ 
rived. 


3.2.1  Results  for  (0)^  Coupons 

for  measurements  of  properties  along  the  fiber  direction, 
we  initially  set  ou.  to  screen  the  (0)^  results  with  Instron 
cons tanc- crosshead  displacement  tests.  By  this  screening  we 
attempted  to  verify  quickly  that  mechanical  response  in  the  fiber 
direction  was  not  sensitive  to  time  effects,  and  thus  reduce  test¬ 
ing  of  these  coupons  to  allow  more  experimental  time  for  testing 
ma r rix- dominated  response . 
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Figure  14  shows  typical  Epp(t)  and  l/^2(t)  results  obtained 
for  the  two  materials  under  RT/wei  conditions.  In  this  case,  no 
provision  was  made  for  maintaining  moisture  in  the  coupons  during 
the  tests,  so  one  might  expect  to  observe  slight  property  changes 
during  the  experimental  time  simply  due  to  moisture  content  changes 
This  was  generally  not  the  case,  as  can  be  noted  in  the  figure-- 
Although  a  fair  degree  of  specimen- to-specimen  variability  was 
noted  for  the  materials  tested  (typically  124-138  GN/m?  for  Epp 
and  0.28-0.33  for  1*12)  ,  there  appears  to  be  no  sensitivity  to  time 
for  the  moisture-conditioned  specimens  tested  at  RT . 

In  Figures  15  and  lb  ,  respectively,  are  shown  the  results 
for  one  of  the  materials,  AS/3502,  when  tested  dry  over  the  en¬ 
tire  range  of  experimental  temperatures.  In  the  case  of  the 
Poisson  ratio  (Figure  16),  there  appears  to  be  a  slight  increase 
with  time  at  some  of  the  temperatures.  This  effect  was  not  re¬ 
producible,  however,  as  can  be  noted  bv  comparing  the  two  differ¬ 
ent  data  sets  shown  for  the  436°K  (325&F)  temperature. 

On  the  basis  of  tests  such  as  these,  then,  we  conclude  that 
the  principal  moduli  Epp  and  ^pg  are  t ime- insens i tive ,  in  agree¬ 
ment  with  other  work  such  as  Ref.  8  for  example.  A  recent  paper 
(Ref.  13)  on  the  time- insensitivity  of  Poisson's  ratio  for  an 
unfilled  adhesive  lends  considerable  support  to  this  viewpoint.. 


3.2.2  Results  for  t-45)?s  Coupons 


tensile  creep  results  on  (d"45)?g  coupons  were  used  to  derive 
compliance,  and  thus  relaxation  modulus  (Ep-(t)),  master  curves. 
From  these  results  for  Ep(t),  shear  modulus  data  were  extracted 
using  the  relationship 


CL2(t; 


2(l  +  j/xv(t)) 


(1) 


Typical  tensile  creep  results  for  one  of  the  materials,  AS/350 
are  given  in  Figure  17.  Using  the  procedure  previously  noted  in 
Figure  9  ,  master  compliance  curves  were  extracted  by  simple  hori¬ 
zontal  shifting  of  the  data. 


We  have  initially  planned  to  provide  an  analytical  descrip¬ 
tion  of  our  data  by  assuming  that  the  compliance  could  be  repre¬ 
sented  by  a  power  law  in  time,  i.e. , 
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D(t)  =  0o  +  Dp  tn 


(2) 


where  D(t)  represents  the  measured  'ensile  compliance,  and  both 
0o  and  Dp  are  temperature  dependent.  Other  previous  work  on 
fibrous  composites  (e.g..  Ref.  7  and  10)  demonstrated  that  such 
a  description  is  quite  useful. 

Repeated  attempts  by  us  to  apply  Eq  .  (2)  howevf,  met  with 
only  limited  success  j.n  that  we  were  not  able  to  find  a  single 
n-value  which  appeared  to  describe  satisfactorily  the  data  at 
all  temperatures.  In  those  cases  where  the  description  appeared 
good,  the  derived  n-value  tended  to  be  quite  small,  e.g. ,  about 
0.05,  whereas  work  on  a  similar  material,  AS/3501-6,  in  Ref.  10 
derived  n^  0.2. 

Lou  and  Schaperv  (Ref.  3)  found  that  recovery  curves  are 
particularly  useful  for  determining  the  exponent  n,  so  we  also 
employed  that  approach.  Due  to  the  rather  weak  t ime- dependence 
of  the  present  materials,  the  recovery  strains  are  quite  small, 
typic’lly  0.001-0.005%  at  RT  dry.  and  there  was  considerable 
scatter  in  our  data.  Although  n  values  around  0.05  were  again 
indicated,  no  consisten.  value  could  be  obtained.  At  the  sugges¬ 
tion  of  Dr.  Knauss  ,  t.hei  we  simply  adopted  the  "visual"  super¬ 
position  procedure  previously  noted.  All  master  curves  were 
thus  formed  by  horizontal  shifting  on  log-leg  plots. 

Results  from  such  a  procedure  for  (9:45) --c  specimens  arc- 
given  in  Figures  18  and  19.  The  data  points  in  Figure  18  are  taken 
from  the  100-minute  dry  creep  tests  for  T300/5208  at  different 
test  temperatures ,  while  the  dashed  curve  represents  similar 
data  for  the  AS/3502  material.  it  is  interesting  to  note  (1) 
the  quite  similar  behavior  for  the  two  materials;  (2)  the  rela¬ 
tively  weak  t ime- dependence  for  dry  materials,  and  (3)  the  near 
linearity  over  roughly  seven  decades  in  time.  The  log-log  slopes 
for  both  materials  are  similar  (about  0.01)  and  in  excellent 
agreement  with  independent  results  obtained  at  Texas  A&M  (Ref. 14) 
on  the  AS/3502  material. 

The  shift  factors  a-p  are  plotted  in  the  inset  to  Figure  18  in 
a  conventional  Arrhenius- type  form,  i.e.,  log  a-p  vs.  reciprocal 
Kelvin  temperature.  If  one  forces  a  straight  line  through  the 
log  a-p  data  to  extract  an  apparent  activation  energy,  the  results 
tend  to  range  about  20-25  k  cal/mole,  which  is  reasonable  for 
these  types  of  materials. 
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This  type  of  straight  graphical  shifting  is  similar  to  that 
used  by  Crossman,  et  al.  (Ref.  15)  on  HMF  330C/934  relaxation 
modulus  data.  If  one  uses  the  roughly  linear  portions  of  Cross¬ 
man's  relaxation  modulus  master  curve,  the  magnitude  of  the  log- 
log  slope  there  is  approximately  0.02,  similar  to  our  value  noted 
above . 

In  order  to  provide  some  assessment  of  these  master  curves 
derived  from  superposition  of  short-term  results,  we  conducted 
two  longer-time  tests  on  the  T300/5208  material.  These  results 
are  shown  as  the  solid  curves  in  Figure  18,  and  represent  a  5000 
minute  test  at  339°K  (150°F)  and  a  7000  minute  test  at  4  Z  2°K 
(300°F) .  These  results  are  positioned  in  time  with  the  log  ay 
values  previously  derived.  The  agreement  between  the  longer 
tests  and  the  superposed  short-time  tests  is  quite  good. 

Turning  now  to  our  results  for  wee  coupons  saturated  at  75T, 
and  95,o  Rh,  we  employed  the  sane  graphical  shifting  procedures 
as  were  used  for  the  dry  data.  We  assumed  here  that  the  moisture 
content  for  a  given  coupon  was  constant  for  all  test  temperatures 
although  there  was  typically  a  moisture  increase  in  the  specimens 
during  the  immersion  tests  of  about  10T  (relative  to  the  initial 
moisture  content).  This  represents  a  rather  small  perturbation 
on  the  overall  results. 

Figure  19  summarizes  the  wet  (and  dry)  results  for  the  T300/ 
5208  material.  The  AS/3502  behavior  was  nearly  identical,  as 
previously  noted.  In  order  to  form  the  combined  uTet/dry  master 
curves,  we  have  assumed  that  moisture-only  and  moisture  plus  tem¬ 
perature  results  are  superposable  in  time  exactly  in  the  same 
manner  as  temperature-onlv  results.  This  implies  that  shifting 
due  to  moisture  ar.d  to  temperature  are  separable  in  the  sense 
that 

Log  ayil  =  log  aT  +  log  aH.  (3) 

Although  there  is  no  a  priori  reason  why  (3^  should  hold, 
the  smoothness  of  our  shifted  curves  and  the  overall  consistancy 
and  reproducibility  of  results  demonstrate  its  usefullness  in 
the  present  case. 

In  Figure  19  we  have  employed  a  reduced  time  scale  which 
extends  roughly  out  to  a  usual  aircraft  service  life,  about  15 
years.  Within  that  time  frame  are  shown  three  different  plots 
of  the  master  curve  for  the  +45  modulus  (reciprocal  of  compliance 
The  top  curve  uses  RT  dry  conditions  as  the  reference  and  demon¬ 
strates  only  small  modulus  changes  over  an  aircraft  lifetime. 
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Realizing  that  no  aircraft  may  be  expected  to  spend  a  service 
life  at  RT  dry  conditions,  the  other  two  solid  curves  in  the 
figure  show  the  shift  in  response  which  occurs  if,  respectively 
380°K  (225°F)  dry  (or  equivalently  RT  wet) ,  or  334°K  (150°F) 
wet  conditions  are  taken  as  reference.  In  the  latter  case, 
we  see  that  rather  dramatic  changes  in  shear  modulus  may  be 
expected  in  times  of  the  order  of  a  month. 

The  inset  to  Figure  19  presents  the  log  aT^  shift  surface 
for  a  range  of  temperature  and  humidity  conditions.  Although 
the  humidity-induced  effects  are  not  dramatic  at  room  tempera¬ 
ture,  the  combined  effects  of  elevated  temperature  and  moisture 
may  be  considerable. 

In  producing  this  shear  modulus  time  dependence,  we  have 
used  Eq .  (1)  and  a  time- insensitive  Poisson  ratio  of  vXy  =  0.72. 
Figure  20  shows  typical  v x„  measurements  from  creep  experiments 
at  various  temperatures.  The  two  separate  measurements  at  350°K 
shown  indicate  that  the  variability  in  results  may  be  at  least 
±3±.  Thus  although  a  slight  upward  trend  in  vXy  may  be  noted  in 
the  figure,  most  measurements  fall  within  or  near  the  standard 
deviation  of  the  mean  there  indicated.  V.'e  have  accordingly  taken 
Pxy  =  0.72  throughout,  and  the  overall  consistency  of  results 
would  appear  to  support  both  thi ;  value  and  its  insensitivity 
to  time.  A  change  in  pXy  of  ±0.1  only  results  in  a  shear  nodulus 
variation  of  about  ±5%,  so  again  the  effect  is  not  great. 

3.2.3  Results  for  (90) ^ 5 

Figures  21  and  22,  respectively,  show  the  E22  master  curves 
derived  from  tensile  creep  compliance  measurements  on  our  two 
material  systems.  The  data  in  either  the  dry  or  the  wet  case 
seem  to  conform  nicely  to  smooth  master  curves,  and  the  shift 
factors  derived  from  graphical  shifting  (horizontal  only)  are 
quite  similar  for  the  two  materials. 

One  interesting  feature  of  the  dry  master  curves  in  Figure 
21  is  their  adherence  to  log-log  linearity  over  the  10-11  decades 
of  reduced  time  represented  there.  The  overall  slope  on  this 
log- log  plot  is  in  the  range  (.008-.  010  !  ,  which  is  quite 
similar  to  the  ±45  master  curve  slope.  Complementary  data  ob¬ 
tained  at  Texas  A&M  agree  in  magnitude  to  within  typically  67, 
and  have  an  equivalent  slope. 
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Figure  22  Wet  E 22  Master  Curves 


If  one  attempts  to  superpose  the  wet  and  dry  data,  however, 
curves  of  the  type  shown  in  Figure  23  result.  One  sees  that 
when  no  vertical  shifting  is  used,  a  severe  difference  both  in 
the  wet  and  dry  master  curves  occurs  at  a  reduced  time  of  about 
10G  minutes.  It  is  somewhat  difficult  to  reconcile  this  result 
with  the  previously  derived  wet/dry  master  curves  for  shear  mod¬ 
ulus  . 


Due  to  the  proclivity  of  90°  cot  ~  :>ns  to  damage,  we  undertook 
an  investigation  to  determine  whether  any  damage  (specifically, 
edge  cracks)  was  evident  in  typical  specimens  of  these  types 
which  had  seen  extensive  histories  of  time  and  temperature  at 
load.  Several  T300/5208  coupons  (both  +45  and  90  types)  were 
edge  polished  and  examined  on  a  bench  metallograph .  Although  all 
the  coupons  had  previously  been  subjected  to  many  100-minute  creep 
cycles  at  various  temperatures,  no  evidence  of  damage  was  found. 

Surface  replicates  of  the  edges  were  then  made  while  the 
coupons  were  subjected  to  small  tensile  loads  in  order  to  open 
any  cracks  which  might  have  healed  during  the  previous  examina¬ 
tion  and  handling.  Again,  no  damage  was  noted.  The  specimens 
were  then  placed  in  a  heated  vacuum  chamber  overnight  to  elimi¬ 
nate  any  residual  moisture  which  might  have  been  absorbed  during 
creep  cycling.  Replicas  were  again  made  under  load,  as  before, 
but  no  cracking  was  evident.  Finally  the  9<~°  specimens  were 
subjected  to  several  creep  cycles  at  stresses  between  20  and  30% 
of  ultimate i,  and  for  times  up  to  16  hours  (at  350°K) .  Periodic 
replicas  failed  to  reveal  any  edge  damage. 

On  the  basis  of  these  limited  damage  investigations,  v?e 
concluded  that  there  were  probably  no  gross  differences  in  the 
specimen  integrity  which  could  reasonably  account  for  the  dif¬ 
ferences  in  the  master  curve  data. 

Dr.  Schapery  of  Texas  A&M  raised  questions  about  possible 
physical  aging  effects  (Ref.  16)  in  these  materials  v’hich  might 
account  for  long  term  stiffness  changes  due  to  free  volume  changes. 
Such  a  mechanism  could  possibly  account  for  the  differences  in 
the  long-term  behavior  of  the  dry  90°  coupons  compared  to  the  dry 
45°,  since  one  expects  the  residual  stresses  in  the  specimen  to 
inhibit  aging  effects. 

We  thus  repeated  the  creep  experiments  on  90°  coupons,  but 
with  one  important  c  iange  compared  to  our  previous  procedure. 

Prior  to  each  creep  run,  the  specimen  was  heated  to  around  450°K 
(350°F),  its  previous  cure  temperature,  with  no  lead  applied. 
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LOG  £  (Minutes) 

Figure  23  Combined  Wet/Dry  E22  Master  Curves 


According  Co  Struck  (Ref.  16),  such  a  procedure  should  effectively 
erase  the  material's  memory  of  its  age,  i.e,,  become  "unaged". 

The  modulus  results  obtained  from  such  specimens  are  shown  in 
Figure  24. 

Two  features  are  very  important  in  Figure  24.  Firstly,  the 
long- time/high- temperature  behavior  shows  a  marked  departure  from 
log-log  linearity,  in  contrast  to  the  dry  data  shown  in  Figures 
21  or  we,  and  in  general  agreement  with  the  ^45  data.  Secondly, 
the  master  curves  for  the  two  materials  shown  in  Figure  24  were 
obtained  with  two  different  sets  of  shift  factors  (log  ap)  .  In 
the  T300/5208  case,  the  master  curve  was  formed  with  G^2-derived 
shift  factors.  For  the  AS/3502,  the  actual  90°  data  shown  were 
visually  shifted  to  obtain  the  best  smooth  curve.  The  discrepancy 
is  evident. 

It  is  not  particularly  satisfying  to  find  that  Gp2(t)  and 
1^22(1)  master  curves  have  different  shift  factors,  if  the  time- 
dependence  is  logically  believed  to  be  matrix  controlled.  Figure 
25  demonstrates,  for  the  T300/520S  material,  that  G]_2~ derived 
shift  factors  are  applied  to  the  E22  results,  then  the  long-time 
results  for  wet  and  dry  (unaged)  coupons  superpose  equivalently 
on  the  same  master  curve  (as  was  the  case  for- wet  and  dry  G}_2 
results) . 

In  Figure  26,  we  repeat  the  T300/5208  data  from  Figure  25, 
and  show  the  individual- temperature  results  for  AS/3502,  shifted 
with  the  same  log  ap  values  as  used  for  the  T300/5208.  It  does 
not  appear  implausible,  given  the  error  bars  indicated  (from  re¬ 
peated  runs)  that  all  the  90°  data  shifts  in  the  same  manner  as 
the  2t45°  results. 

The  apparent  effects  noted  m  these  90°  data  which  we  have 
ascribed  to  physical  aging  certainly  bear  further  investigation. 
Residual  stresses  will  always  be  present  to  different  extents  in 
different  layups,  and  these  stresses  will  be  affected  by  the 
time- temperature-moisture-environment .  It  may  be  very  important 
to  ensure  that  the  effective  ages  of  specimens  used  in  comparison 
tests  are  in  fact  comparable. 

3.2.4  Summary  of  Shift  Factors 

In  Figure  27  are  shown  the  t irne- temperature-humidity  shift 
factors  derived  from  graphical  log-log  shifting  of  the  ^45  results 
and  the  90°  results,  independently.  All  shift  factors  given  there 
are  referenced  to  RT  dry  conditions.  In  each  case  shown,  the 
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Figure  25  Long-Term  Test  Results  Shifted  with  aj  from  M5's 


results  for  the  two  materials  tested  have  all  been  included  to 
obtain  the  fitted  curves  shown  of  the  form 

T 

log  aT11  =  -  -  +  B.  (4) 

We  include  all  six  cur .^s  to  document  our  overall  results. 

As  noted  in  the  previous  subsection,  however,  we  would  tend  to 
recommend  that  the  +45°  coupon  shift  factors  be  used  for  all  data. 
Due  to  the  scatter  of  some  of  the  90°  results,  this  practice  would 
not  be  inconsistent  with  our  overall  results,  and  it  does  offer 
the  more  analy tically-reasonable  result  that  ai1  time  dependent 
results,  being  t-aceable  to  matrix  effects,  are  described  by  a 
single  set  of  shift  factors . 


3.3  Laminate  Properties 

In  our  tests  on  laminate  tensile  compliance,  very  weak  tire 
dependence  was  observed  for  coupons  with  any  fibers  parallel  to 
the  load  direction.  For  example,  in  an  aS/  3502.  (0/±45)5  laminate 
tested  dry,  the  total  change  in  measured  compliance  at  436°K.  was 
about  4%  over  a  100-ruinute  creep  cycle  while,  for  comparison,  a 
similar  test  on  a  t 45  coupon  yielded  a  roughly  12%  change  over 
the  same  time  interval. 

Wet  tests  on  the  fiber-dominated  laminates  yielded  about  the 
same  time- dependent  results  as  the  dry.  In  repeated  runs  under 
the  same  environment  for  a  (0/145) s  laminate,  we  found  a  varia¬ 
tion  of  4%  for  P.T  dry  and  about  5%  for  RT  wet.  Thus,  the  mea¬ 
sured  changes  in  compliance  during  100  minutes  were  typically  the 
same  as  the  test-to-test  variability. 

Tn  Figure  28,  for  example,  are  shewn  dry  test  results  on  a 
(0/90/145) g  laminate.  Little,  if  any  time  dependence  is  noted, 
and  the  change  in  short  time  compliance  from  RT  to  464°K  (375°F) 
is  only  6%.  Although  the  compliance  in  Figure  28  changes  in  a 
regular  way  from  the  lowest  to  highest  test  temperature,  that 
particular  result  was  probably  fortuitous,  since  the  run- to-run 
variation,  as  previously  noted,  was  4-5%.  The  6%  change  indicated 
in  the  figure  is  in  fact  roughly  the  same  as  the  modulus  variation 
noted  in  0°  lamina  tension  tests  over  this  temperature  range  dur¬ 
ing  lamina  design  data  testing  at  General  Dynamics. 

In  order  to  remain  cognizant  of  the  experimental  variability 
and  its  relationship  to  apparent  time  dependence,  in  Figure  29 
are  plotted  'master"  curves  for  two  fiber  dominated  laminates 
under  dry  conditions.  The  data  envelopes  shown  there  include 
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Figure  28  SX(0  for  (C/90/+45)s  T300/5208 


nearly  all  tests  on  these  laminates  at  all  temperatures,  and  the 
elevated  temperature  results  have  been  horizontally  shifted  in 
time  using  the  log  a^  values  derived  from  the  G]_2  results.  Also 
included  in  that  figure  are  predictions  using  measured  principal 
properties  in  an  elastic  lamination  code,  with  the  curves  simply 
calculated  stepwise  in  time.  The  predictions  fall  reasonably 
within  the  scatter  band  of  data  over  the  entire  range  of  reduced 
time  indicated. 

The  experimentally- derived  master  curve,  and  a  similar  pre¬ 
diction,  for  a  matrix-dominated  laminate,  both  wet  and  dry,  are 
shown  in  Figure  30.  Again,  the  shift  factors  derived  from  the 
145  data  were  used.  We  observe  that  the  experimental  master  curve 
is  quite  smooth,  and  the  prediction  is  very  good  over  the  entire 
range.  We  will  discuss  this  particular  laminate  prediction  in 
somewhat  more  detail  in  Section  V. 
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Figure  30  Measured  and  Predicted  Laminate  Response 


3.4  Residual  Strengths 


Residual  strengths  were  measured  under  ambient  conditions  on 
approximately  50  specimens  which  had  been  creep-tested  under  var¬ 
ious  conditions,  and  these  results  are  shown  in  Tables  IV  and  V 
for  the  two  materials . 

There  is  really  little  generalization  which  can  be  made  of 
these  results.  Realizing  that  in  all  cases  the  changes  between 
residual  and  baseline  properties  are  only  slight,  if  any,  the 
comments  indicated  in  Table  VI  might  be  made. 


TABLE  VI  RESIDUAL  STRENGTH  TRENDS 
uayup/'Condition  Trend(Residual  Properties  Compared  to  Baseline) 


(90)  wet 

Some  reduction 

(T45)  dry 

No  change 

(90/ .45 )  dry 

Slight  increase 

wet 

Slight  decrease 

(0/^45)  dry,  wet 

Slight  increase 

(0/90/ ±"4 5)  dry 

Slight  decrease 

wet 

Slight  increase 
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3.5  Stress-Strain  vs.  Strain  Rate 


As  a  slight  addition  to  the  originally-planned  test  pro¬ 
gram,  we  conducted  a  brief  experimental  analysis  of  the  effects 
of  strain  rate  on  tensile  stress-strain  testing.  Instron  stress- 
strain  to  failure  tests  were  run  on  (-45)2c  coupons  of  T300/5208. 
Specimens  were  tested  at  RT ,  380°K  (225°Fj',J  422°K  (300°F),  and 
crosshead  rates  of  0.00085  crn/sec  (0.02  in/min)  ,  .002  (.05), 

.004  (.1),  .0085  (.2),  .02  (.5),  .04  (1),  and  0.085  (2).  Although 

only  one  coupon  per  condition  was  tested,  it  was  felt  that  useful 
trends  might  be  observed,  since  the  inherent  time-dependence  of 
the  material  should  be  expected  to  modify  the  observed  tensile 
behavior  at  various  strain  rates,  particularly  at  elevated  tem¬ 
perature  . 

As  expected,  considerable  scatter  in  the  results  was  evident 
due  to  our  use  of  only  one  coupon  per  condition.  The  following 
trends  in  the  results  are  suggested,  however. 

o  There  appeared  to  be  a  slight  increase  in  stress  ultimate 
as  the  strain  rate  increases. 

o  Strain  ultimates  were  quite  erratic,  but  a  peak  appeared 
around  .004  cm/sec  crosshead  rate,  followed  by  a  valley 
at  somewhat  higher  rates,  then  a  slight  increase  at  the 
highest  rates. 

o  The  Poisson  ratio  appeared  relatively  constant,  al¬ 

though  with  considerable  scatter,  under  all  test  condi- 
t ions . 

o  At  low  stresses,  say  about  13.8  MN/in^  (2000  psi)  ,  the 
longitudinal  modulus  tended  to  increase  as  the  strain 
rate  decreased. 

Although  these  results  were  quite  sketchy  due  to  the  limited 
amount  of  data,  one  should  observe  that  caution  need  be  exercised 
in  the  application  of  tensile  stress-strain  data  for  matrix-dom¬ 
inated  coupons.  Changes  in  measured  properties  do  occur  as  the 
testing  rates  change.  Thus,  for  example,  if  baseline  tensile 
properties  are  to  be  used  to  set  parameters  for  a  fatigue  test, 
the  baseline  rates  and  environment  should  duplicate  the  planned 
fatigue  parameters  as  closely  as  possible. 
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SECTION  IV 


4.  EXPANSIONAL  STRAINS 

In  this  section,  we  present  our  experimental  results  on  the 
longitudinal  and  transverse  strains  observed  in  the  two  material 
systems  due  to  moisture  and  temperature  effects. 


4.1  Moisture  Absorption 

In  order  to  obtain  mois ture- induced  swelling  strains,  uni¬ 
directional  (0)g  specimens  were  conditioned  at  four  constant 
humidities  at  room  temperature.  For  comparison  purposes,  however, 
additional  absorption  measurements  were  conducted  on  these  speci¬ 
mens  at  322°K  (120°F)  and  339°K  <150°F)  and  at  RT  on  samples  of 
all  layups  used  in  this  program.  Samples  for  all  these  uptake¬ 
swelling  measurements  were  typically  3  cm  (1.2  in)  square. 

Specimens  were  conditioned  in  appropriate  temperature  cham¬ 
bers  over  saturated  salt  solutions,  where  the  relative  humidities 
achieved  are  as  indicated  in  Table  VII  below  for  exposures  at 
339°K  (150°F) . 


TABLE  VII  CONDITIONING  PARAMETERS 


Solution 

RH(%) 

KF 

21 

Cr03 

46 

NaCl 

74 

Pb (N03) 2 

93 

The  variation  in  humidity  with  temperature  for  these  solu¬ 
tions  is  not  particularly  significant.  The  lead  nitrate  solu¬ 
tion  possibly  shows  the  greatest  variation,  rising  from  93%  at 
339  K  to  9 7 . 5%  RH  at  RT .  We  will  ordinarily  refer  to  this  high¬ 
est  humidity  used  as  95%,  which  represents  the  average  of  this 
observed  range. 

Moisture  uptake  results  for  these  graphite-epoxy  materials 
are  shown  in  Figures  31  -  33  .  After  initial  drying  to  an 

equilibrium  weight  value,  the  specimens  were  periodically  weighed 
with  a  Mettler  Model  BA28/BL20/CT10  automatic  digital  balance. 
Three  replicates  were  used  for  each  data  point.  For  one  condi¬ 
tioning  case  (RT/9?%  RH)  we  exposed  identical  unidirectional 
specimens  of  both  materials  in  two  different  chambers  to  assess 
conditioning  variabilities.  The  maximum  difference  in  weight 
gain  observed  in  these  "identical"  experiments  was  slightly  over 
10%  (e.g.,  1.0  t  0.1%  moisture  by  weight)  which  is  not  atypical 
for  such  experiments,  even  though  the  standard  deviation  in  resin 
content  is  only  about;  1“ 3% . 


We  may  observe  from  these  figures  that  the  moisture  absorp¬ 
tion  characteristics  for  the  two  material  systems  are  quite  simi¬ 
lar,  although  the  T300/5208  specimens  generally  tend  to  exhibit 
a  slightly  higher  (~ 10%)  weight  gain  than  the  AS/ 3502.  We  also 
note  from  Figure  31  that  the  addition  of  elevated  temperature 
(339°K)  causes  roughly  a  10%  increase  in  moisture  pickup,  even 
though  the  effective  humidity  is  somewhat  lower  at  that  higher 
temperature . 

Figure  34  compares  the  moisture  uptake  under  identical  con¬ 
ditions  for  two  different  6-ply  layups  of  AS/3502.  The  satura¬ 
tion  levels  are  the  same,  but,  as  might  be  expected,  slight  dif¬ 
ferences  in  rate  occur  for  the  differing  stacking  sequences. 

A  summary  of  the  apparent  maximum  moisture  content,  Mm(%) , 
as  a  function  of  relative  humidity  is  given  in  Figure  35  .  We 
have  in  the  past  experienced  difficult  in  obtaining  reliable  and 
reproducible  results  at  very  low  humidities.  The  21%,  RH  data 
shown  there  are  no  expection.  Those  specimens  were  not  weighed 
until  about  250  days  into  the  exposure  in  order  to  avoid  the  long 
re-equilibration  times  noted  after  opening  the  low-humidity  cham¬ 
bers.  Those  results  were  will  inconsistent  with  the  other  mea¬ 
surements,  and  were  not  included  in  the  curve  fit  shown. 
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exposure  times,  which  has  been  noted  by  some  (e.g.  ,  Ref.  IV  and 
18  ),  does  not  seem  apparent  in  our  present  results. 
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T300/5208 


Induced  Weight  Gain  for  (0)g  Coupons 


(DAYS)’/, 

Figure  34  Comparison  ol  Moisture  Uptake  In 
Different  6-Ply  Layups 
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4.2  Moisture- Induced  Swelling 

At  the  same  time  each  of  the  coupon  weight  gains  was  moni¬ 
tored,  the  specimen  dimensions  were  measured  using  4-place  micro¬ 
meters.  This  technique  relies  on  the  use  of  a  fixture  to  hold 
both  the  micrometer  and  the  specimens  in  a  fixed  relationship  to 
one  another,  so  that  reproducible  measurements  are  obtained  on 
the  specimen  at  the  same  reference  position  each  time. 


Swelling  measurements  were  made  on  all  three  specimen  dimen¬ 
sions.  Within  the  limits  of  the  measurements,  however,  no  repeat- 
able  swelling  was  observed  through  the  thickness  (due  to  the  6  ply 
=  .08  cm  dimension).  Along  the  fiber  direction,  a  contraction  of 
typically  0.01  -0.03%  was  usually  observed,  but  such  is  orobably 
not  outside  the  data  scatter. 

Figure  36  shows  the  measured  swelling  in  width  (transverse 
to  the  fiber  direction)  for  both  materials  at  room  temperature 
and  at  elevated  temperature  (322  and  339  K) .  The  RT  data  are 
all  quite  similar  for  both  materials,  and  are  conveniently  des¬ 
cribed  by  the  equation  given  there.  At  elevated  temperature  we 
observe  a  trend  toward  somewhat  greater  expansional  strain  for 
a  given  moisture  content. 

As  is  used  in  such  results,  there  is  no  observable  swelling 
until  a  threshold  moisture  content  is  reached.  The  threshold 
here  observed  is  at  M(%)-  0.1,  although  at  these  very  low  swell¬ 
ing  strains  the  threshold  is  somewhat  indefinite  and  might  be 
assigned  somewhat  higher.  Crossman  (Ref.  9  and  15)  found  the 
threshold  at  0.4%  for  GY70/CE339  laminates,  e.g. 

During  the  course  of  the  program,  we  also  attempted  to  mea¬ 
sure  these  swell i n<  strains  by  two  other  methods,  dens itometric 
and  via  our  thermo- nechanica 1  analyzer  (TMA) .  Since  the  TMA 
utilizes  a  highly  precise  probe  and  LYDT  combination,  v:e  felt 
that  this  latter  method  held  great  promise.  We  were,  however, 
not  able  to  reproduce  measurements  due  to  the  difficulty  of 
repeatable  positioning,  of  the  specimens.  Both  methods  yielded 
unsatisfactory  reproducibility  and  were  abandoned  in  favur  of 
the  micrometer  measurements  described. 


a. 3  Thermal  Expansion 

■'he i  ml  expansion  meausrements  on  unidirectional  (90)  ^5  com¬ 
posites  wi-ie  accomplished  using  a  loaded  column  expansion  method 
(Ref.  19)  with  a  Perk in- Elmer  TMS-1  thermomechanical  analyzer 
(TMA) .  A  schematic  diagram  of  this  test  is  shown  in  Figure  37  - 
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Figure  36  Moisture-Induced  Swelling 


Figure  37.  TMA  Schematic  lor  Loaded- Column 
Thermal  Expansion  Tests 


Samples  o£  dimension  typically  1.2/  cm  x  .64  cm.  x  .20  cm  were 
loaded  parallel  to  the  lcng  dimension  by  placing  20-100  gm 
weights  on  top  of  the  vertical  quartz  probe  u’hich  rested  on 
the  test  specimen.  Changes  in  specimen  length  were  then  mea¬ 
sured  using  a  linearly  variable  differential  transformer  through 
which  the  probe  moved.  The  probe  used  was  a  f i at -bot toned  quartz 
expansion  probe  with  a  diameter  of  approximately  0.25  cm.  Further 
details  of  the  techn ique  and  a  specific  discuss  ion  regarding  mea¬ 
surements  on  moisture-saturated  specimens  may  be  found  in  Ref.  19. 

For  measurements  cf  the  thermal  expansion  coefficient 

parallel  to  the  fibers,  only  one  load/rate  combination  was  used 
in  tne  IMA  apparatus,  namely  100  gm  and  20°K/min.  Figures  38 
and  39  summarize  these  measurements  for  dry  and  wet  coupons  .  V/e 
note  from  those  results  that,  as  expected,  ot  i i  is  negative  for 
both  materials  due  to  the  negative  thermal  expansion  coefficient 
oi  graphite.  The  apparently  large  data  scatter  is  to  be  expected 
however  at  these  high  TMA  sensitivities  (0.1/ie). 

The  solid  curves  in  each  case  represent  a  least  squares  fit 
to  the  data  over  the  range  indicated.  It  is  interesting  to  ob¬ 
serve  the  apparent  slight  differences  in  thermal  characteristics 
of  the  T300  and  AS  fihers:  Of  ^  ^  for  T300  is  about  2-3  times  greater 
for  1300  than  for  AS,  i.e.  ,  typically  -0.6/ve/°K  compared  to 
-0.25 ne/  K.  The  fits  would  also  suggest  that  oi-q  decreases  slight¬ 
ly  with  increasing  temperature  for  T300,  while  the  opposite  behav¬ 
ior  occurs  for  AS  fibers.  Both  curves  suggest  a  dry  glass  tran¬ 
sition  temperature  (Tg)  in  the  480-500°K  range,  which  is  somewhat 
higher  than  that  inferred  from  the  ensuing  figures  for  0*22  ■ 

In  order  to  examine  possible  load  and  rate  effects  during 
our  measurements  of  0(22  (thermal  expansion  coefficient  trans¬ 
verse  to  the  fibers),  we  employed  the  load/rate  combinations  in¬ 
dicated  in  Table  VIII. 


TABLE  VIII  LOAD/ RATE  PARAMETERS 
FOR  a22  MEASUREMENTS 


Scan  Rate  (°K/min)  = 

2.5 

5 

10 

20 

Load  =  10  gm 

X 

X 

X 

X 

=  100  gm 

X 

X 
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OO  -  AS/3502  a  (T)^.00072T  -0.484  I/Xc/CK) 
•  -  T300/5208«  ii  (T)  cr  -0  00047T  -  0.418 
TMA  RUN  AT  100  gm  LOAD  AND  2DCK/MIN 


J _ i _ i _ I _ 1 _ L 

220  260  300  340  380  420 

TEMPERATURE  (CK) 


460  500 


Figure  39  Longitudinal  Thermal  Expansion  Coefficients  a  »  for  Dry  Coupons 
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Figure  40  shows  the  measured  an  (dry)  f°r  the  two  mater- 
ials .  Although  data  for  only  four  load/rate  combinations  are 
shown,  the  other  results  closely  parallel  those  presented.  The 
solid  lines  in  the  figures  are  again  least  squares  fits  to  the 
data  up  to  the  cure/postcure  temperature.  The  c<22  values  are 
quite  similar  for  both  materials  (which  is  to  be  expected  since 
the  resin  systems  are  similar)  and  the  indicated  Tg  for  each 
material  is  in  the  range  of  450-480°K,  at  or  slightly  above  the 
cure  temperature. 

Thermal  expansion  coefficients  ex'll  f°r  wet  materials  were 
run  in  the  TMA  at  20°K/min  and  100  gm  load.  Such  parameters 
have  been  found  to  yield  valid  measurements  on  moisture- satu¬ 
rated  specimens,  even  though  some  drying  occurs  at  the  specimen 
surface  (Ref.  19).  Results  for  the  two  materials  are  shown  in 
Figures  41  and  42  for  specimens  conditioned  at  75%  and  957=  RH, 
respectively.  Two  specimens  of  each  material  were  run  in  the 
TMA.  We  have  plotted  results  for  only  one  experiment  in  each 
case.  The  wet  T  values  at  7 5%  RH  conditions  would  appear  to 
be  in  the  375-385°K  range  (215-233°F) ;  at  957 >  RH,  the  values 
are  somewhat  lower,  namely  360-370°K  (188-206°F) . 

As  noted  in  the  previous  section,  our  test  plan  for  time- 
dependent  measurements  on  wet  specimens  was  based  on  an  expected 
wet  Tg  near  the  boiling  point  of  water.  Thus  our  highest  test 
temperature  for  wet  (creep)  specimens  of  372°K  (210°F)  was  al¬ 
ways  at  or  slightly  above  the  respective  transition  temperatures. 

The  glass  transition  values  obtained  from  Otg2  measure¬ 
ments  are  summarized  in  Figure  43  .  We  have  there  indicated  the 
dry  Tg  to  be  at  the  postcure  temperature,  although  our  results 
may  suggest  a  value  somewhat  higher. 


4.4  Summary  of  Constants 

In  Table  IX,  following,  are  summarized  the  hygrothermal  con¬ 
stants  measured  in  this  task  for  T300/5208  and  AS/3502.  Only  one 
set  of  thermal  expansion  coefficients  are  given  for  wet  (95%  RH) 
coupons,  since  the  results  are  quite  similar  for  the  two  humid¬ 
ities  investigated.  The  mois ture- induced  swelling  parallel  to 
fibers,  /?]_]_,  is  taken  as  zero  from  our  results. 
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TABLE  LX.  HYGROTHERMAL  CONSTANTS 


Coefficient 

Environment 

T300/5208 

AS/3502 

yS22  (|i«/°K) 

K  (%,)  =  0 . 5 

1700 

=  1.0 

2440 

=  1 .4 

2800 

D  (cm  /sec) 

RT 

17.0 

12.0 

«11(>ue/0K) 

Dry  ,  300°K 

t 

o 

ax 

-0.27 

400°K 

-0.61 

-0 . 20 

au(^€/°K) 

9  5  Vo  RH 

-1 

.0 

022  ( °K) 

Dry,  300°K 

23.0 

21  .8 

400°K 

26.7 

25.9 

450°K 

28.6 

28.0 

a22(^/°K) 

95%  RH,  300°K 

28.2 

26.4 

370°K 

34.0 

30.8 

SECTION  V 


5.  ANALYTICAL  METHODOLOGY 

Using  principal  moduli  derived  from  our  compliance  curves 
and  quasielastic  conversion,  we  have  performed  predictions  of 
lamina  and  laminate  behavior.  Since  we  have  shown  that  for  many 
cases  of  practical  interest  the  time- dependence  is  not  particu¬ 
larly  strong,  we  have  stressed  the  use  of  simple  predictive 
techniques  already  familiar  in  composite  elastic  analysis. 

As  an  example,  Figure  44  shows  Halpin-Tsai  predictions  for 
Em,  the  effective  matrix  modulus,  and  for  E22 »  based  on  our  ex¬ 
perimentally  derived  master  curve  (wet  and  dry).  We  have 

used  fiber  properties  essentially  the  same  as  those  given  in 
Ref.  20,  and  summarized  below: 


GfLT 

=3.5  MSI 

(24.1  GN/m2) 

kf-p-p 

=  2.3  MSI 

(15.8  GN/m2) 

Vf 

=  .68 

vt 

=  .20 

^m 

=  .37 

* '  s 

=  1,2 

The  data  points  (o)  in  Figure  41  are  our  experimental  data, 
while  the  solid  curves  are  Halpin-Tsai  results.  The  vertical 
bars  on  the  E22  data  represent  typical  scatter  from  our  90°  cou¬ 
pon  results.  The  E22  data  are  those  obtained  as  described  in 
the  previous  section,  i.e.,  by  using  the  G-j^  shift  factors.  The 
data  points  (□)  are  results  from  Ref. 21  on  heat  3502  resin,  while 
the  (*;  points  are  derived  from  results  in  Ref. 22  for  5208  resin. 

We  have  based  the  Halpin-Tsai  calculations  on  our  G-^  data. 
Given  that  (wet/dry)  master  curve,  and  using  reasonable  values 
for  fiber  properties  and  geometrical  factors,  with  no  unusual 
adjustment  of  parameters,  we  obtain  a  very  reasonable  descrip¬ 
tion  of  E.,?  results.  Additionally,  the  predicted  matrix  modulus 
agrees  wett  with  independent  measurements  on  the  cured  resin 
material . 
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Using  the  time  dependence  tor  and  Co->  shown  in  Figure 

44,  we  have  used  our  elastic  code  to  calculate  the  average  stiff¬ 
ness,  stepwise  in  time,  or  a  wet/drv  (90/^45)s  lav  mate.  lor 
these  calculations  we  used  li p p  -  131  Gfi/m^  (19  MSI),  and  v  ]  •)  " 
0.31,  which  are  "average"  results  from  our  measurements  and  as¬ 
sumed  insensitive  to  time. 


CM 

E 


tSi 

=3 


a 
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Figure  44  Measured  and  Predicted  Lamina  Properties 


Figure  30  in  the  previous  section  presented  the  prediction 
obtained  on  thin  basis,  along  with  the  experimental  master  curve 
for  this  T300/5208  laminate.  The  agreement  between  measurement 
and  prediction  is  within  about  470  over  the  entire  range  of  reduced 
time.  This  would  seem  to  lend  considerable  weight  to  our  conclu¬ 
sions  that  time- dependent  rests  on  t45  coupons,  along  with  Halpin- 
Tsai  predictions  and  conventional  analysis  procedures  provide  a 
good  assessment  of  the  expected  material  response  to  environment. 

Also  presented  previously  in  Section  III  (Figure  29)  were 
predictions  for  fiber- dominated  laminate  response.  In  that  case, 
the  measured  values  for  all  principal  properties  were  used  in 
the  elastic  code  (racher  than  F-22  derived  from  Halpir.-Tsai  com¬ 
putations)  .  Equally  good  predictive  description?  of  the  experi¬ 
mentally-measured  response  were  obtained. 
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For  the  cases  we  have  studied,  i.e.,  for  constant  loads  and 
various  constant  tmepe r a ture/ humidity  conditions ,  it  is  clear 
that  commonly-used  elastic  techniques,  when  used  with  time-varying 
principal  properties,  provide  excellent  descriptions  of  laminate 
response . 

In  order  to  provide  a  more  general  framework  for  analyses , 
however,  we  have  obtained  exponential  series  fits  to  our  time- 
dependent  principal  relaxation  moduli  (as  in  the  case  of  Ref.  9, 
e.g.).  Texas  A&M  University  provided  these  collocation  fits  to 
our  master  curve  results,  and  typical  results  for  the  T300/5208 
(±45)  coupons  are  shown  in  Figures  4  5  and  46. 

In  one  case  (Figure  45)  ,  %-rdecade  spacing  was  used  and  the 
fit  is  excellent.  The  slight  oscillation  at  long  times  could  be 
pushed  off  to  longer  times  by  including  a  "guessed"  extrapolation 
of  the  data.  One  extra  data  point  was  in  fact  used  at  long  time 
by  assuming  the  log-log  slope  is  constant  between  the  actual  last 
two  points  and  an  extra  point  %  decade  beyond  the  data.  The  same 
extrapolation  method  was  used  for  1-decade  spacing  in  Figure  45, 
with  an  extra  point  one  decade  beyond  the  actual  data. 

The  exponential  series  may  be  conveniently  employed  in  heredi¬ 
tary  integrals  for  predicting  laminate  behavior  under  time-varying 
loads  and  environments.  A  detailed  discussion  of  this  procedure, 
along  with  a  general  description  of  viscoelastic  methodology  for 
laminates,  is  presented  in  Appendix  C  ,  "Notes  on  Viscoelastic 
Analysis  of  Laminates",  prepared  for  this  program  by  R.  A. 

Schapery. 
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SECTION  VI 


6.  EFFECTS  OF  NON-UNIFORM  MOISTURE 
DISTRIBUTION 

In  this  task  we  set  out  to  determine  the  effects  of  non- 
uniform  moisture  distribution,  both  experimentally  and  analy¬ 
tically,  with  the  framework  of  methodology  previously  discussed. 
Our  test  coupons  in  this  case  were  24-ply  layups,  (-45)^s  and 
(90/±45>45,  of  both  materials. 

Coupons  were  moisture  conditioned  at  333°K  (140°F)  for 
several  days,  at  which  time  their  total  moisture  content  reached 
about  0.3-0.47,.  Strain  gages  were  then  applied  (and  baked  to 
cure  the  gage  adhesive) ,  and  the  specimens  further  moisture- 
soaked.  Three  coupons  of  this  set  were  then  stress-strain  tested 
to  failure  at  either  RT  or  339°K  (150dF)  for  baseline  values, 
while  another  set  of  three  coupons  were  immersed  in  the  wet 
creep  chamber  at  the  corresponding  temperature.  One-hundred 
minute  creep  tests  were  then  performed  once  per  day  for  1-3 
days.  These  creep  specimens  were  subsequently  tested  (stress- 
strain  to  failure)  for  residual  properties  at  RT.  A  moisture 
traveler  followed  the  coupons  at  all  times  so  that  daily  weight 
measurements  could  be  compared  with  the  predictions  (both  total 
moisture  content  and  moisture  distribution)  of  our  finite  dif¬ 
ference  moisture-diffusion  code. 

Figure  47  demonstrates  the  measured  (data  points)  and  predic¬ 
ted  moisture  content  in  a  set  of  (-45)0S  T300/5208  coupons  during 
this  conditioning  and  test  sequence.  All  specimens  tested  gave 
roughly  equivalent  results  in  terms  of  both  the  experimental 
and  analytically-predicted  data.  We  will,  therefore,  confine 
our  remaining  discussion  to  this  particular  sat  of  coupons  as 
representative  of  all  our  non-uniform  moisture  distribution  data. 

Figure  4b  illustrates  the  analysis  procedure  used  for  these 
data.  Beginning  in  the  upper  left  coraer  and  proceeding  clock¬ 
wise  through  the  four  items  presented,  we  first  observe  the  pre¬ 
dicted  moisture  distribution  through  the  thickness.  Given  this 
M(%)  per  ply,  we  proceed  to  our  relationship  between  moisture 
content,  M,  and  relative  humidity,  <£  ,  determined  in  Section  IV. 
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having  thus  a  relative  humidity  <t>  corresponding  to  each  ply,  we 
can  proceed  to  the  shift  factors,  log  a-j- ,  presented  in  the  lower 
right  as  a  function  of  <j) .  For  the  4>  representative  of  each  ply, 
we  read  off  (from  the  339°K  curve)  the  corresponding  shift  fac¬ 
tor  for  that  ply.  Given  the  log  a-p  values,  w!e  sirrplv  enter  the 
shear  modulus  master  curve  at  the  appropriate  reduced  time,  ob¬ 
taining  a  Gp2  for  that  ply.  This  procedure  is  thus  repeated  for 
discrete  time  steps  during  the  duration  of  the  (creep)  loading. 

Given  then  the  Gg2(t)  appropriate  to  each  ply,  we  utilize 
the  halpin-Tsai  relations  to  compute  Eo2(t),  and  our  previously 
measured  values  for  the  t ime- i nsensi t i ve  moduli  Kjp  and  fp?- 
These  principal  properties  are  then  used  in  a  ply-by-ply  elastic 
lamination  code  to  predict  the  average  tensile  stiffness  as  a 
function  of  time. 

Results  of  this  procedure  are  given  in  Figure  49.  We  plot 
there  the  measured  results  for  modulus  (compl  iance"  •*■)  on  two 
successive  days  (essentially  no  change),  the  previously  measured 
results  for  (.45)  saturated  at  75%  RH ,  and  the  predicted  response 
obtained  as  described  above. 

We  note  that  the  time-dependence  is  similar  in  each  case, 
and  the  prediction  is  within  about  6-7%  in  magnitude  of  the 
measured  value. 

Tab  presents  the  strength  tests  obtained  on  these  cou¬ 

pons  wit'n  non-uniform  moisture  distribution.  The  75%  RH  uniform 
moisture  data  arc  included  for  comparison.  It  is  not  clear  that 
any  substantial  comments  may  be  made  about  those  results. 
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LOG  t  (MINUTES) 

Figure  49  NonUrtiftirm  Moisture  Tests/Predictions 


TABLE  X  KT  STRENGTH  TESTS 


MATERIAL 
&  LAYUP 

MOISTURE 

OONDTTICN 

STRESS 

ULTIMATE 

STRAIN 

ULTIMATE 

% 

WSESfSSMm 

T300/5208 

DRY  BASELINE 

161 

23 

I 

(±45) 

75%  UNIFORM 

160 

23 

HI 

(4) 

(1) 

NON  UNIFORM 

177 

26 

HI 

(1) 

(.1) 

1 

NON  UNIFORM 

— i 

26 

___(*) 

(RESIDUAL) 

| 

(.3) 

— 

AS/ 3502 

DRY  BASELINE 

165 

24 

1.33 

(±45) 

(2) 

(.3) 

(.06) 

75%  UNIFORM 

154 

22 

1.96 

-  -  - 

(3) 

(.4) 

(.47) 

NON  UNIFORM 

171 

25 

2.03 

(1) 

(.2) 

(.12) 

NON  UNIFORM 

174 

mm 

2.07 

(RESIDUAL) 

(2) 

(-39) 

T300/5208 

DRY  BASELINE 

148 

— 

.63 

(90/145) 

(6) 

wm 

(-02) 

75%  UNIFORM 

1 

34 

■na 

H 

(16) 

■n 

NON  UNIFORM 

193 

28 

.81 

(9) 

(1) 

(.05) 

NON-UNIFORM 

32 

.99 

(RESIDUAL) 

| 

(1) 

(.09) 

AS/ 3502 

DRY  BASELINE 

188 

27 

.82 

(90/±45) 

(13) 

(2) 

(.09) 

75%  UNIFORM 

31 

.99 

■I 

CD 

(.07) 

NON  UNIFORM 

239> 

35 

(3) 

(.4) 

Hna 

NGN  UNIFORM 

234 

34 

1.16 

(RESIDUAL) 

(3) 

(.4) 

(.07) 

*  Strain  Gage  Failure 
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SECTION  VII 


7.  CURE  CYCLE  OPTIMIZATION 

The  intent  of  this  task  was  to  use  existing  and  developed 
analytical  tools  in  conjunction  with  our  measured  material  prop¬ 
erties  to  develop  a  procedure  for  minimizing  residual  stresses 
induced  during  the  cure  cycle. 

This  task  was  performed  for  this  program  by  Y.  Weitsman  of 
Texas  A&M  University,  and  a  complete  description  of  his  procedure 
is  given  as  Appendix  D  to  this  report,  entitled  "Minimization  of 
Residual  Thermal  Stresses  in  Cross  Ply  Laminates."  That  work 
presents  an  analytical-numerical  scheme  to  minimize  the  residual 
thermal  stresses  due  to  cool-down  of  a  cross  ply  laminate. 

As  an  example  of  the  application  of  the  procedure  described 
in  Appendix  D,  Dr.  Weitsman  was  provided  with  power  law  fits  to 
our  90  compliance  results  as  described  in  Figure  50.  In  that 
figure  we  show  our  previous  results  for  wet  and  dry  90°  coupons, 
in  which  the  master  curves  for  the  two  conditions  were  obtained 
by  graphical  (horizontal)  shifting  of  the  90°  data.  We  used 
these  results,  rather  than  the  more  consistent  data  set  which 
results  from  using  G^-derived  shift  factors,  in  order  to  examine 
the  widest  possible  range  of  our  data. 

For  each  case  shown  in  the  figure,  the  data  were  fit  with 
the  power  law  form 

D(0  -  V«+Wq  (5) 

and  the  four  sets  of  fit  parameters  given  correspond  to  the  four 
indicated  solid  curves.  Based  on  these  sets  of  parameters,  the 
optimized  cooling  paths  (for  AS/ 3502)  derived  are  shown  in  Figure 
51. 


In  each  case,  the  cool-down  path  corresponds  to  an  initial 
rapid  temperature  drop  from  Tc  (cure  temperature)  to  the  indicated 
Tt_.  follov’ed  by  a  (nearly-linear)  cool-down  to  a  temperature  T2, 
the  cooi-down  being  constrained  to  occur  with  100  minutes  in  this 
case.  All  cases  shown  yield  similar  residual  stresses. 
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For  comparison  purposes,  Table  XI  gives  the  calculated  vis 
coelastic  and  elastic  residual  stresses  for  each  of  the  four 
cases.  The  moduli  used  for  the  elastic  computations  were  given 
bv 


TABLE  XI  RESIDUAL  STRESSES  (MN/nt2) 


Case  No. 


Viscoelastic 


Elastic 


a.  Difference 


11 

.90 

11 

.91 

11 

.99 

12 

.  59 

3  . 19 
3.S6 
.4 .25 


SECTION  VIII 


8.  CONCLUSIONS  AND  RECOMMENDATIONS 

We  have  obtained  a  rather  complete  linearly  viscoleastic 
characterization  of  two  graphite-epoxy  material  systems  by 
measuring  the  principal  properties  as  a  function  of  time,  tem¬ 
perature,  and  humidity.  We  have  used  this  experimental  charac¬ 
terization,  along  with  slightly-modified,  but  conventional, 
analysis  techniques  to  predict  time- dependent  laminate  response. 

We  believe  the  results  are  of  general  applicability  for  ascer¬ 
taining  the  magnitude  of  time- temperature-humidity  effects  on 
structural  composites.  Specific  conclusions  and  recommendations 
are  summarized  below. 

o  Immersion  testing  is  a  viable  procedure  for  short-term 
experiments  on  "wet"  specimens.  The  problem  of  specimen 
dryout  is  circumvented  by  such  procedures,  and  a  reason¬ 
able  choice  of  fluids  for  the  immersion,  particularly 
those  with  high  boiling  points,  allows  simulation  of  any 
relative  humidity. 

o  The  principal  properties  E-^  and  1^2  are  not  sensitive 
to  time,  temperature,  and  moisture.  This  feature  also 
appears  to  hold  for  the  Poisson  ratio  vKy  of  (+45) c 
coupons . 

o  Long-term  results  are  predictable  from  short-term  tests 
by  simple  horizontal  superposition  on  a  log-log  scale. 
Master  curves  extending  over  10  decades  in  reduced  time 
were  produced  by  such  superposition  of  2  decade  tests  in 
different  environments.  Verification  of  such  master  curves 
by  longer  term,  elevated  temperature  tests  extending  for 
nearly  4  decades  in  time  was  accomplished  for  (1:45) c 
results . 

o  Moisture  and  temperature  effects  appear  independently 
superposable  from  the  (+45)  ^  results ,  that  is,  the  com¬ 
bined  shift  factors  may  be  represented  by  a-j-^j  =  a-papj . 

This  is  not  necessarily  a  general  result,  but  in  the 
present  case  appears  to  provide  a  good  overall  descrip¬ 
tion  of  the  data. 
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o  Prediction  of  laminate  master  curves  is  possible  using 
conventional  analysis  procedures.  That  is,  for  the 
loads  and  environments  tested,  quasi-elastic  analysis 
using  stepped  intervals  in  time  to  follow  the  measured 
principal  properties  provides  excellent  predictions  of 
laminate  time-dependent  behavior. 

o  Time- dependent  tests  on  1"45  coupons  for  G-j^Ct)  ,  with  the 
Halpin-Tsai  equations,  provide  an  excellent  assessment 
of  the  expected  material  response  to  environment.  The 
derived  G^2(t)  values,  along  with  reasonable  fiber  and 
geometric  factors,  yield  good  agreement  with  both  our 
measured  E22(ti)  results  and  independent  neat  resin  mea¬ 
surements.  These  G-j^Ct)  and  E 22 (t)»  along  with  measured 
Eh  and  1/^2  values  then  provide  very  acceptable  laminate 
predictions  when  used  with  quasi-elastic  lamination  pro¬ 
cedures  . 

o  Considerable  difficulty  was  encountered  in  obtaining  con¬ 
sistent  results  with  90°  coupons.  Although  damage  effects 
were  suspected,  no  evidence  of  such  was  found.  These 
difficulties  could  perhaps  be  avoided  by  relying  on  (*45)c 
measurements  to  obtain  all  time-dependent  results,  as 
summari  the  previous  item. 

o  Evidence  was  found  for  possible  physical  aging  phenomena 
in  the  90°  results.  Specimens  which  had  "aged"  for  about 
one  year  tended  to  give  different  time-dependent  results 
than  new  specimens.  The  aged  specimens  exhibited  less 
time  dependence  than  expected  on  the  basis  of  +45  results. 
If  the  specimens  were  ’’annealed"  near  the  cure  temperature 
briefly  before  test,  however,  their  creep  response  was 
similar  to  that  observed  for  the  ±45  coupons .  This  is  a 
possible  example  of  the  aging  effects  noted  in  Ref.  16 , 
by  which  the  moduli  of  amorphous  polymers  changes  due  to 
a  reduction  in  segmental  mobility  (decrease  in  free-volume) 
with  time.  If  this  effect  is  indeed  significant,  it 
implies  that  the  usual  superposition  procedures  will  not 
apply.  We,  therefore,  suggest  some  degree  of  caution  in 
measurement  and  prediction  of ^time- dependent  properties. 

One  should  at  least  ensure  that  comparisons  are  made  be¬ 
tween  materials  similarly  aged.  Even  better,  however, 
the  aging  effects  bear  scrutiny  in  some  detail.  Since 
the  materials  are  expected  to  become  "stiffer"  because  of 
such  aging,  procedures  such  as  we  have  here  employed  would 
tend  to  present  a  conservative  prediction  of  long-term 
behavior . 
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Shift  factors  derived  for  (dry)  elevated  temperature 
tests  are  cc rrelatable  with  thermal  expansion  of  the 
materials.  This  suggests  that  the  shift  factors  (dry) 
are  related  to  free  volume.  The  equivalent  correlation 
for  moisture  effects  (swelling  and  shift  factors)  does 
not  seem  as  simple. 

o  All  master  curves  were  derived  by  simple  horizontal 
shifting  on  log-log  scales.  Although  smooth  and  re¬ 
peatable  curves  were  derived,  this  does  not  address  the 
possible  requirements  for  vertical  shifting.  It  is 
clear  that  combinations  of  vertical  and  horizontal  shift¬ 
ing  could  probably  produce  equally  good  master  curves, 
but  the  arbitrariness  of  such  a  procedure  is  question¬ 
able  if  the  data  set  does  not  appear  to  require  it. 

o  Although  limited  tests  were  performed  on  the  rate  de¬ 
pendence  of  tensile  ultimate  properties,  it  is  clear 
that  cognizance  of  such  effects  must  be  maintained.  The 
rate  used  for  tensile  tests  should  be  similar  to  any  rate 
pertinent  to  use  of  the  tensile  properties,  e.g.,  fatigue 
testing . 

o  In  general,  the  time- dependence  of  these  graphite-epoxv 
materials  is  not  particularly  strong,  especially  for  a 
reference  RT-dry  condition,  where  significant  effects 
occur  only  after  many  decades  in  time. 


o  For  other  conditions  mere  representative  of  service, 

however,  such  as  warm/wet  environments,  dramatic  changes 
I  in  matrix-dominated  properties  may  be  expected  in  times 

!  of  the  order  of  months  or  even  weeks.  At  the  verv  least, 

I  i  it  would  appear  appropriate  for  the  designer  to  have 

_  i  access  to  GpgCO  master  curves  for  various  reference  con- 

ditions ,  so  that  at  least  a  preliminary  assessment  of 
j  ‘  potential  property  reduction  might  be  made. 

j  i  o  The  work  reported  here  has  been  concentrated  in  the  re- 

;  ■  -'j  gion  of  linear  viscoelastic  response.  Further  work  should 

\  '  V  be  conducted  to  examine  the  behavior  in  regions  of  non- 

v  linearity. 

;  / 

o  A  general  procedure  for  laminate  analysis  under  condi- 
|  ^  tions  of  transient  and  non-uniform  temperature  and  mois- 

j  *  ture  has  been  presented.  (Appendix  C) 


<  n 


o  Certain  minimization  of  residual  stresses  due  to  cure 
may  result  by  employing  an  optimized  cool-down  path 
after  cure  (Appendix  D) .  Although  the  reduction  in 
residual  stresses  '-re  not  dramatic,  they  serve  to  indi¬ 
cate  that  additional  benefits  may  accrue  if  a  more 
scientific  approach  is  taken  toward  the  entire  produc¬ 
tion  cure  process.  It  is  hoped  that  future  programs 
will  address  this  area  in  more  detail. 
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General  Dynamics  Specif ications  FMS-2023A  establishes  the 
minimum  requirements  for  graphite  -  fiber-reinforced  epoxy  lami¬ 
nating  materials  for  primary  aircraft  structure.  Prepreg  in¬ 
spection  and  mechanical  property  requirements  on  the  supplier 
and  on  GD  for  acceptance  will  be  controlled  by  that  specifica¬ 
tion  and  FPS-2003,  "Advanced  Composite  Materials  Test  Specifi¬ 
cation  .  " 


General  Dynamics  Specification  FPS-2021C  establishes  the 
requirements  for  lay-up,  inspection,  and  cure  of  laminates  using 
FMS-2023  materials.  Ply  inspection,  lay  up,  preparation  for 
cure,  storage,  cure  cycle  temperature  and  pressure  profile  con¬ 
trol,  and  documentation  will  be  controlled  by  that  specification 
and  FPS-2002B,  "Process  and  Quality  Control  of  Adhesive  Bonded 
Boron  and  Graphite  Composite  Assemblies." 

All  laminates  produced  for  this  program  will  be  monitored 
as  outlined  below  to  ensure  the  integrity  of  fabricated  test 
specimens . 

1.  Incoming  Prepreg  Inspection  and  Basic  Material  Properties 

Per  FMS-2023A,  to  include: 

o  Cured  ply  thickness 
o  Tack 

o  Mechanical  properties: 

-  0°  flexural  strength  at  RT  and  450°K  (350°F) 

-  0  flexural  modulus  at  RT 

-  90°  flexural  strength  at  RT  and  450°K 

-  Horizontal  shear  strength  at  RT  and  450  K 

2.  Cure  Cycle  Parameter  Control 
Per  FPS-2021C,  to  include: 

o  Ply  inspection,  layup  and  cure  preparation 
o  Cure  cycle  pressure  and  temperature  profile  control 

Additionally  to  include. 

o  Dielectric  monitoring  of  cure  for  detection  of  irregular¬ 
ities  by  comparison  with  standard  recordings  of  verified 
panels . 
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3.  Cured  Panel  Inspection 


o  Ultrasonic  C-scan  of  panels  prior  to  machining  of  speci¬ 
mens  . 

o  Mark  reject  areas  to  recover  acceptable  areas 
o  Resin  content  measurement  of  ques.ionable  C-scan  areas 

4.  Test  Tabs 

o  Standard  unidirectional  panel  for  mec.anical  ,  roperties 
as  in  (1)  above. 

o  Test  tab  of  each  layup  for  use  in  Task  3,  moisture  trav¬ 
elers  for  Task  1  and  2,  and  property  rerheck,  as  necessary. 

5 .  Pretes t  NDI 

o  Visual  and  optical  examination  of  each  test  specimen  for 

-  General  appearance : 

Scratch  or  scribe  marks,  cracks 
Peel  plv  removal 
Disfiguration  or  distortion 
Unusual  discoloration 
Uniform  thickness  and  width. 

-  Edge  conditions: 

Smooth  and  parallel 

Nicks,  crack,  scratches 

Fibers  displaced  or  pulled  out. 

-  End  tabs 
Proper  taper 

Placement  relative  to  specimen  and  each  other 
Proper  bonding. 

o  Proper  layup  and  number  of  plies 

o  Random  sampling  of  ultrasonic  attenuation  (low-level 
porosity) 

o  Enhanced  x-ray  of  questionable  edge  conditions 

6.  Dimensional  Characterization 

Width  at  three  locations  in  gage  section 
o  Thickn  ss  at  six  locations. 


7.  Microphotographs 

o  60  x  for  voids,  cracks,  delaminations ,  ply  count  and 
orientation  on  at  least  one  specimen/group . 


8.  Fiber  Volume 

o  Specific  gravity  of  cured  laminates 
o  Fiber  and  resin  content. 

9.  History 

o  Complete  material,  panel,  specimen  history  of  above 
measurements . 
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Director,  Mechanics  and  Materials  Research  Center 
Texas  A&M  University,  College  Station,  TX  77840 
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1.  EXPONENTIAL  SERIES  FIT  TO  RELAXATION  MODULI 


Let  D(%)  denote  a  creep  compliance,  such  as  sxx(^)  °r  S22(£»)  ■ 
The  relaxation  modulus  corresponding  to  D(£)  is  denoted  by  EQj=)  ; 
it  can  be  calculated  from  the  following  simple  relationship  with 
negligible  error  (Ref.  1): 


E  <£)  = 


3i  n  mr 


J _ 

DC£) 


(1) 


where 


n  =  d%p 


(2) 


For  all  compliances  obtained  in  this  study,  n  is  sufficiently 
small  to  permit  use  of 

E<0  -  nh)  (3) 

Stress  analysis  of  laminates  is  facilitated  if  all  time-de¬ 
pendent  prinicpal  lamina  moduli  are  expressed  in  terms  of  a  sum 
of  exponentials  in  reduced  time  .  Namely,  let  Ep<  denote  the 
so-called  "Prony"  series, 


E.„fc)  =  E e  +  EE/e  ^ 


(4) 


The  constants  Ee,  E^,  and'C^  can  be  found  by  using  the  collocation 
method  given  in  Ref.  2.  The  quantity  Ee  is  equal  to  the  equili¬ 
brium  modulus , 

Ee  =  E(oo)  (5) 
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which  can  be  taken  as  zero  for  the  range  of  reduced  time  covered 
by  our  data.  The  series  (4)  is  then  equated  to  the  given  (mea¬ 
sured)  relaxation  modulus  E  (%)  at  a  set  of  N  values  of%,  denoted 
by  .  The  time  cons  tants  '"Cp  are  given  as 

~Cl~2Xl,  l  (6) 

Table  C- 1  lists  the  values  of  the  coefficients  and  time  con¬ 
stants  obtained  from  fitting  (4)  to  S22~^($p  for  AS/3502  at  points 
one-half  decade  apart  in  reduced  time.  For  example,  write 

Then  K 2 -j_  are  the  coef ficienizs  listed  in  Table  C-l  .  Evaluation  of 
the  series  at  points  intermediate  to  the  fitted  data  points  serves 
to  check  the  accuracy  of  the  series  representation;  we  have  found 
excellent  agreement  with  the  given  function,  such  as  S2?~'L(£) 
which  is  plotted  in  Figure  C-l. 

Before  turning  to  laminate  stress  analysis,  it  should  be  noted 
that  the  exponential  series  (4)  enables  one  to  calculate  easily 
the  complex  modulus,  which  is  useful  in  vibration  analysis. 
Specifically,  let  E"  denote  the  complex  modulus 

~E+jE  (8) 

where  j  =  V^T.  Then  corresponding  to  an  exponential  series  fit 
to  the  relaxation  modulus, 
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TABLE  C-l  Exponential  Fit  Parameters  For  E22  Master  Curve  (AS/3502) 
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TABLE  C-l  (Continued) 

Exponential  Fit  Parameters  For  E22  Master  Curve  (AS/3502) 
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(9) 
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uu'-T'El 

ZJrSi 


and 

y  do) 

c  Ji1 

where  w  is  the  frequency  and  the  constants  (Eg.  E^,  tL)  are 
those  found  previously  from  the  reciprocal  of  creep  data. 


2.  STRESS  ANALYSIS  OF  VISCOELASTIC  LAMINATES 
2.1  Lamina  Constitutive  Equations 
By  employing  the  accurate  quasi-elastic  method  described 
in  Ref.  1,  the  relaxation  reduced  lamina  stiffness  and  the 

laminate  stiffness  can  be  derived.  Following  the  notation  in 
Ref.  3,  the  so-called  reduced  stiffnesses  are 


0  = 

(U) 

S^s- 

o  -  ~  S'^ 

c  c  _  e4 

(12) 

0  =  — ^ — 

^zz  s  s  -s2, 

°II  az 

(13) 

°66 

(14) 

where  S^j  =  S^ j (£)  are  the  principal  creep  compliances  for  the 
lamina.  Exponential  series  can  be  fit  to  the  =  Q  —  (£)  by 

the  collocation  method  described  above. 


The  full  set  of  lamina  constitutive  equations  in  the  lamina 


i  4 


principal  coordinate  systems  is 


-  Q.fc-Ot.AT  6,Am}  + 

"^12.  =  Qee  JV,*} 


(15) 


where  Qj_j  f  }  is  shorthand  notation  for  the  hereditary  inte¬ 


gral 


(16) 


and  "I"  represents  any  one  of  the  arguments  (e.g.,  £. ,  -oc  ,AT  ,-^AM) 


appearing  in  (15).  Also 
•k 


*.-/£  ,  u/ 


t' 


dtl 

Qtm 


0  ^th  0-  (17) 
where  aj;.j  =  axM(T,M)  is  the  temperature-  and  moisture-dependent 
shift  factor.  For  time- independent  moisture  and  temperature, 


4-  £  , 


(1.8) 


The  expansion  coef  f  icients  CX_i ,  0C2  ,  <Sp  ,  and  ySg  >  in  general 
may  depend  on  reduced  time  (Ref.  1).  If  thev  are  independent  of 
reduced  time,  but  are  functions  of  temperature  and  moisture, 
they  are  properly  interpreted  as  secant  (not  tangent)  values . 


10? 


Also,  for  this  case,  it  can  be  shown  that,  without  any  loss  of 
generality  oc^  can  be  taken  as  the  coefficients  in  a  reference 
moisture  state;  i.e.,odj_  =ot^(T),  but  in  general  ./Si  -/SiCT.H). 

2.2  Laminate  Analysis  for  a  Spacewise  Uniform, 
Transient  Temperature 

For  the  case  in  which  the  temperature  and  moisture  are  either 
time- independent ,  or  spacevise  uniform  or  both,  elastic  solutions 
by  means  of  the  quasi-elastic  method  can  be  used  to  predict  lami¬ 
nated  stresses  and  deformations  (Ref.  1).  Each  ply  is  assumed 
to  depend  on  the  same  shift  factor  a^. 

Let  us  illustrate  the  method  by  showing  how  laminate  response 
can  be  calculated  for  the  case  of  a  t ime- dependent ,  spacewise- 
uniform  temperature,  in  which  stresses  result  from  thermal  ex¬ 
pansion  (or  contraction)  and  external  mechanical  loading. 

As  the  first  step,  we  solve  the  problem  for  an  elastic  lami¬ 
nate  using,  for  example,  the  equations  in  Ref.  3;  in  this  step, 
the  mechanical  loads  and  temperature  change  AT  are  to  be  time- 
wise  constant.  However,  the  reduced  stiffnesses  for  each  ply 
are  to  be  those  in  Eqs .  (11)  -  (14),  which  depend  on  reduced 
time.  This  so-called  quasi-elastic  solution  should  be  evaluated 
at  a  series  of  reduced  times,  one  or  one-half  decades  spacing  of 
the  reduced  time  will  normally  suffice  if  an  interpolation  is 
used  when  performing  the  integrations  described  below. 


Suppose  the  inputs  are  the  forces  Mx,  Nv,  '?xy  and  moments 
Kx,  My,  and  Hxy  •  and  temperature  change  ZiT.  Consider  any  one 
response  quantity  of  interest;  say  the  stress  (<5"2)k  in  the  k-th 
ply.  The  quasi-elastic  solution  ;nn  always  be  written  in  the 
form, 

(<£.)*  =  A.IV  X5f\+  AfeM7  *  A7at 

(19) 

where,  in  general,  A  ^  are  functions  o  reduced  time,  they  are 
simply  the  coefficients  found  from  elastic  analysis  using  time- 
dependent  properties  m  Eqs.  (11)  -  (14).  For  time-dependent 
inputs , 

(ci)k=  ^°) 

where,  as  before,  Ail}  is  shorthand  for  the  integral 

xL[z}= 

0  (21) 
Thus  the  problem  is  reduced  to  elastic  analysis  plus  numer¬ 
ical  or  analytical  evaluation  of  the  integrals  (21).  A  numeric- 
ally-efficient  method  for  evaluating  such  integrals  whenX^  is 


expressed  in  an  exponential  series  is  described  below.  For  iso¬ 
thermal  problems,  the  integration  can  often  be  accomplished 
analytically,  depending,  of  course,  on  the  time-dependence  of  the 


2.3  Analysis  for  Transient:  and  Nonuniform  Temperature 

and/or  Moisture 

For  the  general  class  of  problems  considered  here,  quasi- 
elastic  analysis  cannot  be  used  beyond  the  calculation  of  the 

Qij  in  Eqs .  (11)  -  CIA) .  Instead,  an  incremental  time-marching 
numerical  method  of  solution  is  required,  such  as  used  in  Ref. 

4  for  cross-plv  laminates . 

A  very  efficient  numerical  solution  for  solving  analogous 
problems  involving  isotropic  viscoelastic  materials  was  devel¬ 
oped  several  years  ago  (Ref.  5)  and  is  detailed  in  Ref.  6.  It 
is  based  on  the  use  of  moduli  which  are  represented  in  exponen¬ 
tial  series  in  reduced  time.  The  shift  factor  may  vary  from 

point  to  point.  For  ease  in  reference,  it  is  reproduced  here, 
with  an  appropriate  change  of  notation  to  treat  anisotropic 
materials.  The  method  is  particularly  suited  to  the  finite  ele¬ 
ment  method.  Basically,  the  lamina  constitutive  equations  (15) 
and  (10)  are  reduced  to  elastic  equations  for  the  current  time 
step. 

In  order  to  illustrate  the  method,  it  will  be  sufficient  to 
consider  the  numerical  evaluation  of  an  integral  which  is  repre¬ 
sentative  of  those  appearing  in  Eq .  (15): 

°  (22) 
Let  <£>--¥  (23) 


112 


where 


0*(£)  =  E  0me 


Then,  Equ.  (22)  becomes 

•L 


S=  pe  +  j0”(K-tn&Jt' 


The  discussion  deals  explicitly  with  the  effect  of  tempera¬ 
ture  in  the  equation.  Moisture  can  be  included  in  an  analogous 
fashion . 

2.3.1  Incremental  Determination  of  the  Constitutive  Law 

Subdivide  the  time  span  £ojtl  into  N  intervals  such  that 

to  =  o,  t]_  «  tD  +At1 . t.N  =  tN_x  +  AtN  =  t.  For  S(x,tN), 

£(x,  ty)  write  S  and  C  respectively . 

Assuming  that  S  (x  o)  =  £(x  o)  =  o,  Eq  .  ''‘'5'  c'n  be  written 

Sn  ■ 


Dif-.dt' 

V,  (26) 

Now  within  each  time  interval  approximate  the  time  rate  of 
change  of  each  strain  component  using 

i  -t  ±  ts  iu 

dt  /  o  7  \ 


J  * 


Thus,  Eq .  (26)  becomes 


SH=(Zb£N  */  0^M-S')||'d- 

°  f*' 

|  0*(£N-S')d£' 


\-£  /I 

•6m  i 


4 


Define  the  quantity JJL$  as 

r^-N 


0*«„-£)dt'  +  ft 


and  let 


U =  at„  ‘  €»-> [  0'*(^-K')<a' 

'"Tw-I 

0*  (£„-£')  H 


or,  using  (29)  in  (30) , 

i  f  C  NJ:*(  ll / 


Using  the  definitions  (29)  and  (36) ,  the  hcr.ditary  integral 


The  analog  of  Eq .  (32)  appears  in  elasticity  when  considering 
a  material  subjected  to  thermal  loading  (L^) .  represents  the 

modulus  for  material  point  x  at  time  tj’  and  reflects  time- tempera¬ 
ture  dependence  of  relaxation  behavior  during  the  time  interval 
tN-l~t~tN'  is  c^e  contribution  of  complete  past  strain  and 

temperature  history  on  current  ( t^)  response.  The  procedure  for 
calculating  the  quantities//^  ,  Lr^  is  discussed  below. 

2.3.2  Determination  of 

Consider  Lq .  (29)  and  substitute  in  Eq .  (29)  for0  : 

U 

.  I  fyn  -i  ~Prr,(ZN'K) 

Atw  d£ +0o 

>'■  UA  (33) 

LW-| 


i 
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1 
!  S 


!'  i 


1 1 


The  0m  are  time- independent  material  parameters;  thus 

0*  Ik  f Vmr£“"£,)  di'j + & 

~tu~[ 


(34) 


Defining 


T  -  -L 

°Nm 


J. 


we  have 


H 

JJ-N  a  0^^/m  +  0C 


(35) 


(36) 


The  problem  of  calculating JjL ^  is  reduced  to  that  of  evaluating 
the  integral  J^m .  The  procedure  for  finding  JNm  is  delayed  at  this 
point,  and  attention  is  now  directed  to  the  term  LN  of  Eq .  (32), 

which  is  also  shown  to  depend  on  the  integral  of  Eq .  (35). 

2.3.3  Evaluation  of 

Defining  $*(£.„-  £')  df '  (37) 

o 

Eq .  (31)  can  be  written  as 

LN  ~  [  —  (/46n  “  0a)  £*j-i  1 

7  (38) 

Substitute  Equ.  (24)  in  Eq ,  (37)  for  0*  and  move 0m  from  under  the 
integral  as  before 

tL  -  *»-i 

M  /  <->V-  11/ 

i  &t'  dt 

(39) 


ZN-  E  0m(  re-M£“'£,)  if-dt' 

(71 » I  t  O 


Let 


-W-I 
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Nm 


-W£n-£')  Be 

er  at'dE 


(40) 
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1  ■-  "TTrawrrr 


Hence 


How  expand  as  the  sum  of  two  integrals 


C 


n  tw-z 


Nm 


e  azdt 


-h 


tu-1 


at 


3 1' 


or 


Cwm=  e*"  V"5**  U  H  dt' 


at 


tn-2 


From  definition  Eq  .  (40  it  is  seen  that 
Uz 

/  -/W a£  jr  =  r 

J  0  s€  dt 

o 

Thus,  Eq .  (43)  can  be  written  as 

_  -)W£n-£mJ) 


£  =  0 

^Nm 


✓Tu-i 


N-l/m 


J 

"tw-2 


-/Syt,(£n_,-£;')  _36_  ,/? 

c  at  J 


(41) 


(42) 


(43) 


(44) 


(45) 


Let  ^  =  A£,^  and  use  the  approximation  Eq  .  (27)  in  Eq  .  (45) 


C  =  0 

Eniti  ^ 


'^m4bu 


p^v-i 

*Y 

t-N-Z 


N-i,vn 


+ 


m] 


(46) 


But  the  integral  appearing  in  Eq .  (46)  is  defined  previously  in 
Eq .  (35)  ,  hence 
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(47) 


To  summarize  this  section,  the  terra  L*j  which  represents  a 
contribution  of  past  s train- temperature-mois ture  history  to  Eq . 
(22)  valid  at  time  t^,  is  found  using  Eqs .  (38),  (41),  and  (47)  . 
From  the  form  of  Eq .  (47),  it  is  seen  that  the  effect  of  history 
is  considered  by  means  of  recursive  computation  based  on  certain 
parameters  known  from  the  previous  solution  time  point.  In 
terms  of  numerical  computation,  this  means  that  total  history 
need  not  be  retained  for  calculating  a  constitutive  law  governing 
current  response. 

In  order  to  complete  the  calculations  for  and  Lft  a  pro¬ 
cedure  for  evaluating  the  integral  (35)  must  be  developed, 

2.3.4  Evaluation  of 


Recalling  Lq .  (35) 
ntu 


Jutr]  At  '  J 


Now  expand  Eq .  (35) 

X,  -  srj + 

C  tu-i 

■4„.,  "  j> 

r 

~r  _  1  C  i 

^  L  iu 

.  f  "  -/SnACt-t)!  1 , 2 


^roCK.-K) 


Let  the  change  in  reduced  time  during  the  time  interval 

At^  be  related  to  At^j  through  the  parameter,  a, 

K-K-,  - 


°r  At 

Define  the  quantity  J^tn  as 

r»^N 


(50) 


X-ii  J "  [e  -  e  *“*■*>}  <k 


(51) 


Thus  Eq .  (49)  becomes 


T  =  1 .  J.  .  r 

^tN  a/5m  L  °  J 


In 

+  T* 

y.  WNm 

t-u-i 


(52) 


Using  Eq.  (48)  and  performing  the  evaluation  at  the  limits  of 
integration,  Eq .  (52)  reduces  to 


T  -  . 

[i-e*"^]  -j* 

4m 

A  At  ^Nm 

P*n 

(53) 

Finally  consider  J^m  defined  in  Eq .  (51) 

X X \  [e'/8"c<"'£)- 

Hi.! 


(54) 


Let 


+ 


At. 


At 


N 


which  is  to  say  that 

A  -X 

~  <ALH  _  , 

7T-— 1 fr-Q+t*;  tH.^t<t„ 


n; 


(55) 


Affecting  a  change  of  variable  in  Eq .  (54) 


T*  -  J_  f  r 
'-'urn  Z  J 


]  dt 


(56) 


e  '  -  er 

Eq.  (56) is  in  a  form  suitable  for  numerical  evaluation  using 
Gaussian  Quadrature.  One  final  step,  however,  must  be  developed 
in  order  to  complete  the  determination  of  Eq .  (56).  Reduced 
time  and  real  time  must  be  functionally  related  over  the  interval 


tN- 1  -  t  -  CN- 

2.3.5  Evaluation  of  £ 

At  a  given  point  in  space  approximate  the  temperature  as  a 
linear  function  of  time-''  within  the  interval  At^; 


*  The  order  or  this  approximation  is  consistent  with  that  assigned 
for  each  strain  component,  see  Eq .  (27). 


Define 


Q(T> 


efCT  ’  dT' 

•* 


(63) 


Thus,  using  Eq .  (63),  Eq .  (62)  becomes 

£■  ^N-I  +  Auj  QCTn-,)J  (64) 

for  £ui 

Noce  that  the  change  in  reduced  time  in  the  interval  Atj]  is 
found  from  Eq .  (64)  to  be 


a4n=a7n[Q(tJ-Q(v,)J 


(65) 


Suppose  that  the  temperature  does  not  change  in  the  interval 
Atjg,  then  reduced  and  real  time  are  proportional.  Thus,  when 
a2N  =  (64)  and  (65)  become 

£  *  Sn-1  +  (64a) 


A^a  eHT""\tN 


(65a) 


2.3.6  Summary  of  Computational  Steps 

The  formulae  presented  in  the  preceding  sections  are  for 
evaluating  the  hereditary  law  at  time  tN  for  a  linear  visco¬ 
elastic  material:. 


Assuming  that  a  solution  has  been  obtained  at  time  t^_^  and 
that  the  temperature  and  moisture  distribution  is  known  at  tjj , 
then  we  can  proceed  to  evaluating  the  response  at  t^  as  follows 

AT* 


^X-l 


TnV^- 

in” 


At 


A 


N 


ZN  At, 
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Recalling  the  approximation  for  T(t)  within  Atjj 


T  -  A]_j^  +  A2^t , 

Then  reduced  time  is  calculated  using 

£  -  K., + [Q(t>  -gkt„.,)] 

[if  Am--0,  then  eflT“"'l0i-tN..')J 

^=£-~.,+Aito-QM 

[if  Ajk  =  0 ,  then  +  AtNeffT“'>J 

Q(J)  -  [  efCT  )  ir 

T* 

The  shift  function  0(t)  =  ef^T^  is  conveniently  defined  in 
a  numerical  table  (0i  vs.  T-;  ).  Then  Q(T)  can  be  calculated  using 
numerical  integration  at  all  points  in  T  given  by  the  input  table 
(T*  is  arbitrarily  assumed  to  be  the  lowest  temperature  specified). 
The  value  Q(T)  for  T  t  Tj_  is  selected  by  interpolation  once  Q(T^) 
a rc  established. 

Recti  1  ling  that 

“  Kh-i 

and  that  the  relaxation  function  in  distortion  is  known  in  terms 
of  a  Prony  Series 

M 

zi*  Edme*" 

rt\  *| 

we  can  evaluate  the  integral 


and 


where 


I 


using  Gaussian  Quadrature 
Next  form 


v-^Nno 


jj  -  e^’"A^w] 


KX 


Kirn 


which  allows 


M 

Mh~  0® 

X  M»l 


to  be  calculated. 

Using  the  response  of  the  previous  solution  time  interval,  form 

C  -e^-rc  +f£  -£  'IT  1 

which  leads  directly  to 

XN  =  f  0„Ck 

,v  m»l 


'Nm 


which  in  turn  allows 


Ln-  [tn  -  (/xN-0o)eN.,] 


to  be  evaluated. 


The  constitutive  law  governing  each  component  of  stress- 
strain  is  defined  at  this  point  for  tpj  .  Thus, 

SN  + 

An  equation  of  this  type  is  to  be  derived  for  each  integral 
in  Lq .  (15).  The  problem  is  thereby  reduced  tc  elastic  analysis 
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Minimization  of  Residual  Thermal  Stresses 
In  Cross  Ply  Laminates 


by 

Y.  Weitsman* 


Abstract 


This  work  presents  an  ana  lytical -numerical  scheme  to  minimize  the 
residual  thermal  stresses  due  to  cool-down  of  cross  ply  laminates. 

The  analysis  considers  the  time-dependent  behavior  of  composites 
and  the  calculations  are  based  upon  recent  data  on  the  thermoviscoelast.ic 
response  of  graphi te/epoxy  laminates.  It  is  shown  that  substantial 
reduction  in  the  residual  stresses  can  be  achieved  by  following  an 
optimal  cool-down  path. 


it 


v- 
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1. 


Introduction 


This  work  strives  to  apply  some  recent  data  on  the  time-temperature 
response  of  composites  [1]*,  [2]  to  achieve  a  reduction  of  the  residual 
thermal  stresses  in  composite  laminates  which  arise  upon  cool-down  from 
cure  temperature  to  room  temperature,  by  controlling  the  cool-down 
process. 

An  analytical-numerical  scheme  is  presented  whereby  an  optimal 
cool-down  path  is  determined  which  minimizes  the  residual  stresses. 

The  analysis  is  based  upon  a  recently  developed  scheme  [3]  -  [5]  and 
the  numerical  computations  are  conducted  for  material  parameters  that 
are  representative  of  realistic  properties. 

The  calculations  are  performed  for  symmetric,  cross-ply  laminates 
which  contain  equal  numbers  of  0°  and  90°  plies.  When  such  laminates 
are  subjected  to  fluctuating  temperatures,  the  stiffer  fibers  within 
each  ply  restrain  the  transverse  deformation  of  the  adjacent  laminas, 
and  these  mutual  geometric  constraints  introduce  substantial  residual 
stresses  in  the  laminate. 

Since  temperature  causes  thermal  expansion,  it  can  be  considered 
as  a  stress  inducing  agent.  On  the  other  hand,  temperature  influences 
most  significantly  some  of  the  material  properties  of  fibrous  laminates 
in  a  way  which  enhances  the  relaxation  of  stresses.  The  optimal  time- 
temperature  path  achieves  the  best  interplay  among  the  above  mentioned 
competing  effects  so  as  to  minimize  the  residual  thermal  stress  at  the 
termination  of  the  cool-down  phase. 

2 .  EJ asti c  Analysis 

Consider  a  symmetric  cross-ply  composite  laminate,  consisting  of 
an  equal  number  of  0°  and  90°  plies,  exposed  to  fluctuating  ambient 
temperature  T(t).  For  typical  fiber- reinforced  composites  it  is 
permissible  to  discard  the  precess  of  heat  conduction  and  assume 
spatially  uniform  temperatures  at  all  times  t.  Let  C^,  and 


^Numbers  in  brackets  indicate  reference  listed  at  the  end  of 
this  work. 
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denote  ply  moduli  and  a^,  coefficients  of  thermal  expansion,  where 
subscripts  L  and  T  refer  to  the  longitudinal  (i.e.  parallel  to  fiber) 
and  transverse  directions,  respectively. 


Consider  first  a  geometrically  constrained  laminate  which  is 

subjected  to  a  uniform  temperature-increment  AT.  Elementary 

considerations  yield  the  following  expressions  for  the  stress  o  and  o  : 

x  y 


0  _  90  _  ,r  ,  r  xaT 

ox  -  oy  -  - (C, a,  +  C, ,aT)AT 


1  L  V'12"T' 
'TaT  +  C12  V 


0  -  90  _  !r  .  r  ,.T 

Oy  ox  “ ( Ct<it  +  rt !  )i.T 


0) 


In  (1)  the  superscripts  "0"  and  "90"  refer  to  the  0°  and  90° 
laminas,  respectively,  assuming  that  the  longitudinal  (i.e.  fiber) 
direction  of  the  0°  ply  coincides  with  the  x  direction  of  the  under¬ 
lying  x,  y  Cartesian  system.  (Fig.  1). 

Neglecting  edge  effects,  the  above  stresses  are  relieved  by 

spatially  uniform  strains  e  and  c  so  as  to  yield  null  resultant 

x  y 

forces  F  and  F  in  the  >:  and  y  directions, 
x  y  J 

Due  to  the  assumed  symmetric  layup  of  equal  numbers  of  0°  and 
90:i  plies  wo  have 


The  strains  (2)  give  rise  to  the  following  stresses: 
0  90 


y 


=  (cL  +  c12k 


(3) 


90  _  0_/r  r  , 

o  -o  ■(C-r  +  0lo)r 

x  y  T  12' 

Combining  (i)  and  (3)  and  rec|uiring  null  resultant  forces  v/e  get: 

(4) 

(C,-.,  *■  CT'VT  +  C19ci,  +  C19r>T)AT  +  (C,  +  CT  +  2C12)c  -  0 


0  ,  90  0  x  90 

O  +0  =0  +0 

xx  y  y 


whereby 


C |  |  a-  CTu-|-  C ,  9  ( u |  ^  ) 


(b) 


CL  +  CT  +  2C12 
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0  _  90  _  90  _  0  _  aT  ~  aL 

O  0  “  0  “  "O  “ ’  f  ,  p  '  "n  «> 

x  y  x  y  ST  +  SL  -  2i>12 


3.  Viscoelastic  Formulation 


(7) 


It  has  been  observed  [1],  [2]  that  the  compliances  and  Sp  remain 
constant  at  all  temperatures  belov/  the  glass  transition  temperature 
T  .  In  addition,  it  was  possible  to  express  the  t i me  dependent 
transverse  compliance  in  a  "power  law"  form  as  follows 


S-  =  ST  +  S  tq  (8) 

ilc 

where  $j  is  the  instantaneous  comp! iance,  9  and  q  arc  material 
constants,  and  t  is  time. 

At  elevated  temperatures  Sj  can  be  expressed  Iv,  a  thoriaorheol ogi ca  1 1 
complex  relation 

St(T)  -  :>,(T)  +  Sc[t/o  (T  )]q  (9) 

In  (9)  a(T)  is  the  temperature  dependent  "shift-factor"  function. 

In  the  sequel  we  shall  employ  the  average  value  of  Sj  over  the 
range  of  temperatures  relevant  to  our  problem,  thus  approximating 
(9)  by  a  "therm.-.. rheological  ly  simple"  expression. 

Consequently,  for  the  viscoelastic  counterpart  of  (7),  we  express 
the  compliances  ( S -j.  a  -  2S^)  in  the  form 

B(t)  -  iij  a  Bctq  (10) 

whore 

CT  -  S,avg-  +  S,  -  2S,  and  B  -  S 
II  L  12  c  c 

The  computation  of  stresses  is  greatly  facilitated  when  a 
relaxation  function  E(t)  is  employed  in  place  of  the  compliance  B(t). 
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These  two  material  properties  are  related  by  reciprocal  relations  in 
the  Laplace  transformed-time  domain  ^ 

PE(p)  =  [pB(p)]"1 


Furthermore,  an  approximate  inversion  gives,  to  high  degree  of 
accuracy,  the  expression  ^ 


E(t) 


1 

pB(p) 


(ID 


Substitution  of  (10)  into  (11)  yields 


E(t)  = 


_ 1 _ 

Bj (1  +  Atq) 


(12) 


where 


A  =  2qr(q  +  1'  -3q 


In  the  subse<|uent  analysis  we  shall  further  approximate  (12)  by 
the  following  expression 

E(t)  -  i:0(t.0  +  t)~q  (13) 

where  Fq  and  tg  are  chosen  in  such  a  manner  that  (13)  provides  good 
agreement  with  (12)  over  the  time-span  of  interest. 

The  advantage  of  form  (13)  is  that  it  is  amenable  to  an  analytic 
solution  for  the  optimal  cool-down  path  which,  furthermore,  is  given 
by  elementary  functions. 

Employment  of  (13)  to  represent  the  denominator  in  (7)  gives  the 
following  expression  for  the  stress  o(t),  accounting  for  transient 
viscoelastic  behavior: 


c(t) 


t 

ds 

a(TTsT) 


dAT  (jJ 
dr 


dr 


In  (14)  o 


(14) 
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Consider  now  the  cose  of  cooling  from  the  initial  cure  temperature 
Tj  down  to  room  temperature  in  a  'given  time  interval  t^.  We  search 
for  an  optimal  temperature-time  path  T(t)  which  minimizes  a(t^). 

It  has  been  shown  [3],  [5]  that  this  optimal  path  consists  of 
abrupt  drops  in  temperature  at  t  =  0  and  at  t  =  t^.  The  initial  drop, 
from  T  to  Tg,  is  given  by  solving  the  the  transcendental  equation 

T  -  T  =  ^0-T  05) 

0  1  aTT0T 

In  the  open  interval  0<t  <  tf  the  optimal  path  T(t)  was  shown  to 
be  smooth  and  continuous,  where  it  is  governed  by  the  nonlinear  integro 
differential  equation 


dT  .  r(Mf)  a  ‘  (T(t) )  (16) 

dt  '  rit.tyy  a (rroi o 'Oto-) 

In  (15)  and  (16)  primes  indicate  derivatives  with  respect  to  the 
argument.  Also 


E(tr  t2)  -  E 


ds 

ol'T  (?)) 


With  given  c(t),  a(T)  and  Tj,  the  value  of  TQ  is  determined  by 
(15)  and  the  optimal  temperature  path  can  be  constructed  by  iteration 
through  the  employment  of  (16). 

The  iteration  scheme  consists  of  dividing  the  time  interval  tf 
into,  say,  n  equal  sub- i nterva 1 s  At,  and  selecting  a  guess  value  for 
T(t*)  denoted  by  Tf^.  Then,  at  t  =  t'f  (16)  yields 

dT  j  _  E" (0)  a'(Tf) 

dt  !,  _  -  "  ‘  E  'To)  a(TfT)av'(T’f )  (17) 


The  discretization  of  tf  into  increments  At  gives 
,  F" (0)  a' \Tf) 

U  f  '  Al;  "  'f  "  E '  (6)  a_(r 


Now  that  both  and  T(t^  -  At)  are  determined  it  is  possible 
to  extrapolate  (16)  back  to  t  =  t^  -  2?X  and  determine  Tvt^  -  2qt). 

The  procedure  extends  backward  up  to  t  =  0.  Obviously,  it  is  to  be 
expected  that  the  value  of  T(0),  thus  extrapolated,  will  not 

correspond  to  the  value  of  Tn  as  determined  by  (15).  If  T(0)  >  Tn 

u  (?)  ill  U 

we  must  select  another  guess  value  ,  and  vice  versa. 

The  sequel  of  guess  values  for  T^  are  then  iterated  among  themselves 

until  a  value  of  T^.  is  found  which  yields  T(0)  =  Tq  to  within  a 

desired  accuracy. 

It  is  clear  that  the  appropriate  value  of  T^.  will  generally  diffo 
from  the  room  temperature  T^.  Consequently,  a  second  abrupt  variation 
in  temperature  is  needed  at  time  t^  to  reach  the  prescribed  value  of  T 

With  the  optimal  temperature-path  thus  determined,  the  stresses 
at  all  times  t,  and  in  particular  at  t  =  t|,  can  be  determined  by  the 
numerical,  time-discretized,  evaluation  of  (14). 

A  particular  simplification  occurs  for  a  "power  law"  response,  as 
expressed  in  (13),  when  the  shift-factor  function  is  given  by 

a(T)  --  e/p(-  l  *  (!)  (18) 

In  this  case  it  has  been  shown  [5]  that  the  optimal  path  is 
given  by 


T  1  -  L 

1  C  K  +"  ~Ct 


where 


>■  •  L  "  ‘  ->pCct0) 

while  K  is  determined  as  the  root  of  the  transcendental  equation 


KeCT0  = 


1  q  +  1 
eRt0(K  +  Ctf)q 


Substitution  of  (18)  and  (19)  into  (14)  yields,  after  several 
manipulations : 


t 

r 

r 

1  I 

1 

o(t)  =  (aL  -  <*t)Eq 

L 

/ 

fl< 

K  +  Ct)  q  -  (K  +  Ct)  q 

+  t0 

1 

0 

L 

■md 

r  -  --  -I' 

- 

-q 

+ 

Q  K  q  -  (K  +  Ct)  q 

+  lo 

<Ti  -  V 

i 

L  — 

- 

(20) 


In  (20)  Q  - 


q±l 

(q  +  1 )Ae’B  Lq 


1 


|  'I 

i 


*i 


At  t  r  t^,  a ( )  is  obtained  by  setting  t^  in  place  of  t  in  (20). 
At  t  =  tt: 


o(t+)  =  o(t^)  +  (c<L  -  aT)E0tQC1(Tf  -  Tr) 


(21) 


5.  Numerical  Computations  and  Results 


The  numerical  evaluation  of  (20)  was  obtained  by  discretizing 
the  time  t.  For  this  purpose  the  interval  tf  was  divided  into  n  equal 
sub-intervals  At  =  t^/n  and  n  intermediate  values  of  time 
tj  =  (j  -  1 ) A t ,  j  =  l,2,...n+l  were  considered. 


The  temperatures  T.  =  T(t)  are  obtained  directly  from  (19). 
J  J 

The  stress  =  o(t^)  =  o(0T),  immediately  after  the  initial 


drop  from  Tj  to  T^,  is  given  by 


(22 ) 


while  the  stresses  at  all  subsequent  times,  during  cool-down,  are  given 
by 
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a 


j+i 


) 


(aL  aT)E0 


/  -  -  -i\ 

j 

[Q  y  (K+Cti+i )  q  -  (K+Ctj+1)  qJ 

i  =  l 

L  K  +  cti+i 

4  [o((^Cti)  q  -  (K4Cti+1)  <*))+  t0]‘q 


K+Ct. 


^Q(K"-  (K*Ctj+, )'  «  )  ♦  t0]'q  (T,  -T0) 


(23) 


with  j  =  1 ,2,. . .n. 


The  stress  at  tp  following  the  cool-down  from  T  =  1  ln  _ — 

•  t  C  K  + 


t.o  Tr,  is  obtained  front  (21) 


Ct, 


Numerical  computations  were  carried  out  for  the  following  material 
parameters : 

(a)  Compliances  and  moduli: 

SL  =  0.5277  x  10‘7  (psi)"1  S1?  =  -0.1847  x  10"7  (psi)'1 


ST  :  V  Sctq  Sj  -  5.634  x  10'7  (psi) 


-1 


q  =  0.18 


Sc  "  0.235  x  10  7  (psi)  ^  per(rnin) 


.  \  0 . 1 8 . 


These  compliances  were  fitted  into  a  relaxation-modulus  given  in 
(13),  with  LQ  -  3.18  x  107  psi  and  tQ  =  5.711  x  107  min. 
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(b)  Coefficients  of  thermal  expansion: 

aL  =  -0.17  x  10"6  in/in/l°K  aT  =  7.94  x  10'6  in/in/l°K 

Tj  =  450°K(=350°F)  TR  =  297°K(=75°F) 

(c)  Shi  ft- factor  function 

a (T)  =  exp{-  ^  +  B)  T  in  degrees  Kelvin,  A  =  6.284 
B  =  47.265. 

Results  are  shown  in  Figs.  2  and  3.  In  Fig.  2  the  optimal 
temperature  paths  T(t)  and  the  ensuing  stresses  a ( t )  are  shown  for 
cool-down  times  =  31  min.  (solid  lines)  and  =  100  min.  (dashed 
lines). 

Note  the  discontinuities  n  the  paths  T(t)  and  o(t)  at  t  ~  0  which 
are  independent  of  t^. 

A  second  discontinuity  occurs  at  the  termination  of  cool-dov/n  when 
both  optimal  paths  sustain  sharp  drops  which  lower  the  temperature  down 
to  Tr  =  297r,K.  These  drops  are  accompanied  by  sudden  increases  in 
u(t).  Observe  also  the  typical  steepening  of  the  path  T(t),  and  the 
accompanying  stress  path  o(t),  for  the  shorter  value  of  tf.  We  also 
observe  that  while  T ( t )  decreases  monotonical ly  during  cool-down  the 
stress  o(t)  oscillates  at  short  times. 

Fig.  3  demonstrates  the  reduction  of  residual  thermal  stress 
o(t^)  with  cooling  time  t^.  Note  that  reductions  of  about  25%  are 
attained. 

6 .  Summary  and  Co nclusions 

The  analysis  presented  herein,  which  draws  upon  relevant  data 
for  graphite-epoxy  laminates,  indicates  that  substantial  reductions 
in  the  residual  thermal  stress  due  to  cool-down  can  be  achieved  by 
following  an  optimal  cooling  path  in  time.  Several  approximations 
were  employed  in  modeling  the  visco  =  elastic  response  of  the  laminate. 
These  include  the  averaging  of  the  initial  compliance,  the  utilization 
of  an  approximate  form  for  the  relaxation  modulus,  and  the  assumption 
that  data  from  aged  composites  applies  to  freshly  cured  laminates.  Those 
approximations  and  assumptions  will  merit  more  careful  consideration  in 
the  future.  , 
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Figure  _T ijjes: 

Fig.  1:  A  Symmetric  Laminate  Consisting  of  Equal  Number  of  0°  and  90° 
PI i os . 

Fig.  2:  Optimal  Cool-Down  Paths  T(t),  and  the  Accompanying  Values  of 

the  Thermal  Stress  c(t).  Vs.  Time  t  During  Cool-Down,  for  Two 
Values  of  Cool-Down  Tine,  tf  =  31.6  min.  (Solid  Lines)  and 
tf  =  100  min.  (Dashed  Lines)  T  in  °K  and  o  in  psi .  Note  the 
Sudden  Drop  at  t  =  0  from  Tj  =  45C°  K  to  Tq  =  443.72°  K 
Accompanied  by  a  Sudden  Pise  in  Stress  to  65.1  psi. 

Fig.  3:  The  Residual  Stress  o(tf),  at  Termination  of  Cool-Down, 

Along  Optimal  Path,  Vs.  Cooling  Time  tf  (in  Minutes). 

Note  the  Logarithmic  Scale  for  tp.  o  in  psi.  The  Elastic 
Value  of  1900  psi  is  shown  by  a  Dashed  Horizontal  Line. 
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